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Abstract 


Direct  sequence  pseudo-noise  encoded  spread  spectrum 
waveforms  are  examined  for  their  vulnerability  to  unautho¬ 
rized  intercept  and  message  decoding.  The  interceptor  is 
assumed  to  be  operating  without  the  knowledge  of  the  spread¬ 
ing  code  used.  No  attempt  is  made  by  the  interceptor  to 
duplicate  the  code  in  the  analysis. 

A  general  waveform  equation  is  introduced  to  be  used 
in  the  analysis.  Two  generic  modulation  formats,  designated 
as  format  A  and  format  B,  and  their  associated  transmitters 
are  presented,  and  it  is  shown  that  format  A  and  B  both  fit 
the  general  waveform  equation. 

An  intercept  receiver  model  is  proposed  that  will  de¬ 
code  transmissions  that  fit  the  form  of  the  general  waveform 
equation.  The  receiver  is  analyzed  for  its  success  in  de¬ 
coding  both  generic  modulation  formats.  It  is  shown  that 
for  format  A,  decoding  of  the  message  without  knowledge  of 
the  pseudo-noise  code  is  not  possible  for  this  receiver  when 
only  one  transmission  is  received.  However,  when  two  simul¬ 
taneous  transmissions  are  received,  message  decoding  is  pos¬ 
sible.  Factors  that  affect  successful  decoding  in  this  case 
are  the  carrier  phase  differences  between  the  two  transmit¬ 
ters,  the  offset  between  the  pseudo-noise  code  bits  at  the 
receiver,  and  additive  noise.  It  is  also  shown  that  for 
modulation  format  B,  decoding  of  the  message  is  possible 
when  only  one  transmission  is  received.  Factors  that  affect 


successful  decoding  in  this  case  are  additive  noise,  Doppler 
phase  shift,  and  the  extent  to  which  the  baseband  modulation 
was  linear. 
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INTERCEPT  VULNERABILITY  OF 


DIRECT  SEQUENCE  PSEUDO-NOISE 
ENCODED  SPREAD  SPECTRUM  WAVEFORMS 

I  Introduction 

Spread  spectrum  (SS)  communications  systems  are  those 
which  are  intentionally  designed  to  transmit  information 
over  a  bandwidth  much  wider  than  would  normally  be  required. 
Since  their  beginnings  in  the  early  1950's,  SS  systems  have 
found  both  commercial  and  military  applications.  Some  of 
the  properties  that  have  made  the  use  of  SS  systems  desirable 
are  as  follows  (Ref.  6 * 6 ) s 

1)  Selective  addressing  capability; 

2)  Code  division  multiplexing  for  multiple  access; 

3)  Low  density  power  spectra  for  signal  hiding; 

4)  Message  screening  from  eavesdroppers; 

5)  High  resolution  ranging; 

6)  Interference  rejection. 

The  interference  rejection  property  makes  the  use  of  SS 
communications  systems  attractive  for  use  in  a  high  density, 
multisignal  environment  where  strong  interfering  signals 
would  render  a  conventional  narrowband  system  useless.  An 
SS  system,  in  theory,  will  demonstrate  a  robust  immunity  to 
strong  interfering  signals  and  jamming. 
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Problem  Statement 


The  designer  of  an  SS  system  intends  that  unauthor¬ 
ized  eavesdropping  be  extremely  difficult.  A  method  used  in 
the  past  for  unauthorized  demodulation  has  been  for  the  in¬ 
terceptor  to  duplicate  the  spreading  code  used,  and  synchro¬ 
nously  demodulate  with  this  code.  It  has  been  shown  in  the 
literature  (Ref.  16  and  17)  that  for  a  linearly  generated 
periodic  code  of  length  2n-l  bits  per  period,  the  code  can 
be  reconstructed  by  observation  of  any  2n  consecutive  bits 
of  the  code.  This  method,  however,  can  be  made  prohibitive¬ 
ly  impractical  by  using  non-linear  codes  (Ref.  17).  The 
purpose  of  this  thesis  is  to  show  that  for  a  class  of  modu¬ 
lation  formats,  intercept  and  demodulation  of  SS  waveforms 
is  possible  without  knowledge  of  the  pseudo-noise  spreading 
code,  under  appropriate  conditions. 

Scope 

This  thesis  will  consider  only  direct  sequence  gener¬ 
ated  SS  systems.  The  level  of  presentation  is  directed  to¬ 
ward  the  working  communication  engineer.  Extensive  mathe¬ 
matical  derivations  are  an  unavoidable  necessity.  Although 
an  intercept  receiver  is  proposed,  no  testing  or  experimen¬ 
tation  will  be  accomplished.  The  intercept  receivers  opera¬ 
tion  is  examined  and  performance  predictions  are  made. 

Assumptions 

The  following  assumptions  were  made  to  narrow  the  prob¬ 
lem  to  one  that  could  be  handled  within  the  time  allotted* 
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1.  The  interceptor  is  dedicated,  knowledgeable,  and 
properly  equipped  for  SS  signals.  He  knows  the  general 
frequency  and  modulation  characteristics  of  the  transmitter, 
but  not  the  pseudo-noise  keying  codes  used.  He  can  adjust 
his  hardware  to  fit  a  particular  situation. 

2.  The  interceptor  is  not  interested  in  just  detect¬ 
ing  the  presence  of  an  SS  signal,  but  rather  he  desires  real¬ 
time  message  reconstruction. 

3.  Linearly  generated  codes  are  not  used  because  of 
their  ease  of  reconstruction  by  the  interceptor. 

4.  The  interceptor  is  limited  by  practical  consider¬ 
ations  to  a  "simple"  receiver,  that  is,  one  that  does  not 
entail  a  great  deal  of  sophistication  such  as  the  use  of  a 
massive  computer  center  for  message  processing. 

5.  The  interceptor  is  situated  in  a  multisignal 
tactical  war  environment. 

Approach  and  Sequence  of  Presentation 

The  approach  used  in  this  thesis  is  to  examine  the 
effect  of  processing  direct  sequence  modulated  SS  waveforms 
through  non-linear  elements:  square  law  and  fouth  law  de¬ 
vices.  Although  there  are  many  ways  to  attack  the  intercept 
and  decoding  problem,  this  particular  approach  was  used  be¬ 
cause  performing  squaring  and  quadrupling  operations  on 
pseudo-noise  code  modulated  waveforms  generates  terms  that 
are  independent  of  the  code.  The  results  of  these  non¬ 
linear  operations  are  then  examined  for  possible  exploita¬ 
tion.  The  sequence  of  presentation  is  as  follows.  First,  a 
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chapter  is  devoted  to  an  introduction  to  the  basic  fundamen¬ 
tals  of  SS  systems.  Next,  a  general  waveform  equation  is 
introduced  that  encompasses  several  different  modulation 
formats.  The  general  waveform  equation  takes  pulse  shaping 
of  the  data  into  account  for  offset  quaternary  phase  shift 
keying,  minimum  shift  keying,  and  sinusoidal  frequency  shift 
keying  modulation  schemes.  Two  generic  spread  spectrum  mod¬ 
ulation  formats,  designated  as  format  A  and  B,  and  their 
associated  transmitters  are  introduced  for  eventual  analysis 
of  their  vulnerability  to  intercept.  The  generic  modulation 
formats  are  shown  to  fit  the  general  waveform  equation. 

In  the  next  chapter,  an  intercept  receiver  model  is 
proposed.  The  receiver  is  analyzed  with  the  general  wave¬ 
form  equation  as  its  input.  The  receiver  is  divided  into 
general  sections  designated  as  the  Carrier  Reconstruction 
Loop,  the  Format  A  Decoder,  and  the  Format  B  Decoder.  The 
sections  are  analyzed  in  that  order.  It  is  shown  that  re¬ 
gardless  of  the  modulation  format  or  pulse  shaping  used,  a 
signal  is  developed  in  the  output  of  the  quadrupling  circuit 
of  the  Carrier  Reconstruction  Loop  that  is  free  of  pseudo¬ 
noise  encoding.  Coherent  carrier  references  are  derived 
from  this  signal,  and  it  is  shown  that  in  the  case  of  modu¬ 
lation  format  B,  the  signal  contains  the  desired  message 
information.  Next,  it  is  shown  that  for  modulation  format  A, 
decoding  the  message  without  knowledge  of  the  pseudo-noise 
code  used  is  not  possible  for  this  receiver  when  only  one 
transmission  is  received.  However,  when  two  simultaneous 


transmissions  are  received,  message  decoding  is  possible 
under  appropriate  conditions.  Successful  decoding  in  this 
case  is  dependent  upon  the  carrier  phase  differences  between 
the  two  transmitters,  the  .arrival  offset  between  the  pseudo¬ 
noise  code  bits  at  the  receiver,  and  additive  noise.  These 
effects  are  included  in  the  output  equations.  Reconstruc¬ 
tion  of  the  message  from  the  output  equations  is  demonstrated 
for  ideal  conditions  (i.e.,  negligible  noise,  synchronized 
carriers,  and  no  code  bit  offset).  Finally,  it  is  shown 
that  for  modulation  format  B,  decoding  of  the  message  is 
possible  even  when  only  one  transmission  is  received.  Addi¬ 
tive  noise,  Doppler  phase  shift,  and  the  linearity  of  the 
baseband  modulation  affect  the  success  of  decoding  this  for¬ 
mat. 
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II  Fundamentals  of  Spread  Spectrum 
Communications  Systems 

The  term  "spread  spectrum"  is  applied  to  any  of  a 
number  of  signaling  techniques  in  which  the  transmitted  wave¬ 
form  occupies  a  bandwidth  far  in  excess  of  the  baseband  sig¬ 
nal  (Ref.  9)  and  where  the  envelope  of  the  power  spectral 
density  of  the  transmitted  waveform  is  not  related  to  the 
data.  These  techniques  generally  fall  into  four  categories 
(Ref.  ?:1): 

1)  Direct  sequence  modulated; 

2)  Frequency  hopping; 

3)  Time  hopping; 

4)  Pulse-FM  or  chirp. 

Direct  sequence  signals  are  generated  by  modulating  a 
carrier  with  some  code  sequence.  In  this  type  of  system,  the 
incoming  information  signal  is  generally  transformed  into  some 
digital  format  and  then,  assuming  binary  data,  it  is  added 
modulo  two  to  a  higher  speed  code  sequence.  The  combined  in¬ 
formation  and  code  are  then  used  to  suppressed-carrier  modu¬ 
late  an  RF  carrier.  The  high  speed  code  sequence  dominates 
the  modulating  sequence  and  determines  the  transmission  band¬ 
width  which  is  "spread"  in  comparison  to  the  information 
bandwidth. 

In  a  frequency  hopping  system,  the  carrier  frequency 
is  moved  around  in  frequency.  The  frequency  used  at  any  par¬ 
ticular  time  is  determined  by  a  code  sequence.  This  code 
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sequence  is  the  same  type  that  is  used  in  direct  sequence 
systems  except  that  the  code  clock  rate  is  generally  much 
smaller  because  of  technological  considerations. 

In  a  time  hopping  system,  the  transmission  time  and 
period  are  controlled  with  a  code  sequence.  Time  hopping 
can  be  viewed  as  pulse  modulation  under  code  sequence  control. 
Often,  time  hopping  is  used  in  conjunction  with  frequency 
hopping. 

Pulse-FM  or  chirp  is  a  technique  commonly  used  in 
radar  systems.  This  technique  does  not  use  coding  to  spread 
the  transmitted  signal  bandwidth  as  the  others  do,  but  uses 
swept-frequency  pulses.  The  receiver  employs  a  dispersive 
filter  to  match  the  swept  signal  and  compress  it  into  a  much 
narrower  time  slot. 

This  thesis  will  be  limited  to  consideration  of  direct 
sequence  modulated  SS  systems  because  of  their  more  common 
usage  at  the  present  level  of  technology. 

Spread  Spectrum  Communications  System  Model 

Figure  1  illustrates  the  basic  components  of  a  direct 
sequence  modulated  SS  communications  system  using  binary 
phase  shift  keyed  (BPSK)  signaling.  The  message  waveform, 
m(t),  is  a  bipolar  binary  sequence  with  m(t)  =  ±1.  (Note 
that  multiplication  of  a  carrier  by  plus  one  or  minus  one  is 
equivalent  to  shifting  the  phase  by  zero  or  rr  radians  respec¬ 
tively.  )  The  phase  modulated  carrier  is  multiplied  by  the 
spreading  waveform,  c(t),  a  high  rate  binary  sequence  common¬ 
ly  referred  to  as  a  pseudo-random  (PR)  or  pseudo-noise  (PN) 
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Figure  1.  Spread  Spectrum  Communications  System  Model. 


sequence.1  These  designations  are  a  consequence  of  the  fact 
that  the  sequence  generates  sample  statistics  similar  to  a 
random,  noiselike  waveform.  However,  the  sequence  c(t)  is 
completely  deterministic  and  periodic.  For  extremely  high 
signaling  rates,  the  spectrum  of  the  sequence  is  so  spread 
and  the  spectral  density  so  low  that  to  the  casual  observer 
it  appears  to  be  background  noise.  The  spreading  waveform, 
c(t),  takes  on  values  of  ±1  in  order  to  accomplish  the  re¬ 
quired  phase  shifting. 

The  receiver  has  a  replica  of  the  spreading  code  which 

is  multiplied  by  the  received  signal  to  strip  the  code  off  of 
2 

the  signal.  The  signal  is  then  processed  thru  a  correlation 
detector  whose  output  is  the  sufficient  statistic,  e.  From 
the  sufficient  statistic,  a  decision  is  made  as  to  which  of 
two  hypothesis  is  true.  Under  one  hypothesis,  HQ,  m(t)  =  1. 
Under  the  other  hypothesis,  H^,  m(t)  =  -1. 

It  is  clear  that  the  difference  between  this  system 
and  conventional  ones  is  multiplication  by  the  PN  sequence, 
c(t),  at  the  transmitter  and  receiver.  To  understand  the 
effect  of  c(t)  on  the  system,  the  characteristics  of  PN 
sequences  must  be  understood. 


1The  rate  of  the  spreading  waveform,  c(t),  is  typi¬ 
cally  on  the  order  of  megabits  per  second,  while  the  rate  of 
the  message,  m(t),  is  typically  on  the  order  of  kilobits  per 
second. 

2Note  that  c2(t)  =  1. 
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Characteristics  of  PN  Sequences 

Consider  a  random  binary  sequence.  Such  a  sequence 
could  be  generated  by  repeatedly  flipping  a  "fair"  coin  and 
recording  the  results,  "one"  for  a  head,  "minus  one"  for  a 
tail.  The  resulting  random  sequence  will  have  the  following 
three  "randomness  properties"  (Ref.  10:10  and  15:279) « 

1)  The  Balance  Property.  The  number  of  ones  will  be 
approximately  equal  to  the  number  of  minus  ones. 

2)  The  Run  Property.  Among  the  runs  of  consecutive 
ones  and  of  zeroes,  about  one-half  of  the  runs  of  each  kind 
are  of  length  one,  one-fourth  of  each  kind  of  length  two, 
one-eighth  of  length  three,  and  so  on. 

3)  The  Correlation  Property.  The  expected  value  of 
the  autocorrelation  function  of  the  sequence  would  be  maximum 
at  the  origin  and  decrease  rapidly  away  from  the  origin.  PN 
sequences  exhibit  all  three  of  these  randomness  properties. 

A  PN  sequence  is  made  up  of  bits  of  duration  T  and, 

C 

to  be  consistent  with  the  model  under  consideration,  will  be 
considered  to  have  an  amplitude  of  ±1.  It  is  assumed  that 
the  sequence  is  continuous,  that  is,  there  is  no  space 
between  bits?  that  the  positive  and  negative  bits  occur  with 
equal  probability;  and  that  successive  pu'.ses  are  statisti¬ 
cally  independent.  Figure  2. a  illustrates  a  portion  of  a 
typical  PN  sequence,  c(t). 

The  autocorrelation  function,  R  (t),  of  the  determin- 

c 

istic  PN  sequence,  c(t),  is  the  time  average,  strictly  over 
all  time,  of  the  sequence  multiplied  by  a  time  shifted 
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Figure  2.  Autocorrelation  Function  and  Power  Spectral 
Density  of  a  PN  Sequence,  (a)  Typical  Code  Sequence. 

(b)  Autocorrelation  Function,  (c)  Power  Spectral  Density 
(Ref.  12). 
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version  of  itself  (Ref.  12).  Although  a  PN  sequence  is  peri¬ 
odic,  it  is  assumed  that  the  period,  T,  is  greater  than  the 
duration  of  the  mission.  Thus,  for  all  practical  purposes, 
the  PN  code  period  can  he  considered  infinite.  Then,  the 
autocorrelation  function  is  given  hy  (Ref.  12) 

i  T 

.  R  (t)  =  lim  ■=*  S  c(t)  c(t+r)  dt  (1) 

c  T-«o  2T  _T 


where  t  represents  the  amount  of  the  shift. 

From  the  randomness  properties  of  PN  sequences,  it  is 
clear  that  for  |t|>T  ,  the  average  value  of  c(t)*c(t+r)  is 
zero  since  its  instantaneous  value  has  an  equal  probability 
of  being  +1  or  -1.  For  I  “T !  — T  ,  the  average  value  over  one 

V  ' 

bit  of  c(t)*c(t+r)  is  given  by  the  area  of  the  overlapping 
part  of  the  displaced  bits,  shown  shaded  in  Figure  2. a,  di¬ 
vided  by  T  .  This  area  is  given  by  T  —  f r I .  Thus 
^  c 


RC(T) 


< 


J  (tc-|tI) 

o 


for  It  I  ^  Tc 
for  It  I  >  T 

C 


(2) 


R  (t)  is  plotted  in  Figure  2.b. 

C 

The  power  spectral  density  (PSD),  S  (co),  for  a  PN 
sequence  is  found  from  the  autocorrelation  function,  R_(t), 

C 

using  the  Wiener-Khintchine  theorem* 


S  (w)  =  /  R.(T)  e~^  dT 

U  -00  ^ 


(3). 
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Since  Rc(r)  is  an  even  function,  this  takes  the  form 

00 

S  (to)  =  X  R  (t)  cos(tor)  dr  (4) 

—oo  ^ 


Then  inserting  Rc(t)  and  exploiting  symmetry  gives 

T 

S  (to)  =  I  X  c( T  -  It ( )  cos(toT)  dr 
c  ic  0  c 

sin2(coT  /2) 

=  t  _ - _ 

°  (COT  /2)2 

v 

S  (to)  is  shown  in  Figure  2.c. 

C 


(5) 


Performance  of  Spread  Spectrum  System 

Referring  to  Figure  1,  the  transmitted  waveform  of 
the  SS  model  is  given  by 

s(t)  =  V2E/T  m(t)  c(t)  cos(tOQt)  0  <  t  £  T  (6) 

where  T,  the  message  bit  interval,  is  some  integer  multiple 
of  Tc,  the  code  bit  interval  (T  =  NTQ).  The  received  signal 
is  given  by 


r(t)  =  s(t)  +  n(t)  (7) 

where  n(t)  is  zero  mean,  white,  Gaussian  noise  with  PSD  of 
Sn(to)  =  Nq/2,  and  is  independent  of  the  code,  c(t).  The 
received  waveform  is  multiplied  by  a  stored  replica  of  the 
code  to  give  the  following  (assuming  perfect  synchroniza¬ 
tion)  i 


r'(t)  =  ->/2E7t  m(t)  c2(t)  cos(coQt)  +  c(t)  n(t) 


=  V 2E/t  m(t)  cos(d)0t)  +  n’(t)  (8) 

The  autocorrelation  function  of  n’(t)  is  given  by  the  fol¬ 
lowing  equation: 

R,(W  =  n’(t2)} 

=  ECcC^)  c(t2)}  ECnCtj^)  n(t2)} 

=  Rc(W  (N(/2)  (9) 

Because  of  the  delta  function,  this  has  value  only  when  t^ 
equals  t2>  Then 

R'Cti-tg)  =  Rc(0)  (Nq/2)  6(t1-t2) 

=  (Nq/2)  6(trt2)  (10) 

Thus,  n'(t)  is  itself  white,  Gaussian  noise.  The  problem  is 
thus  reduced  to  determining  whether  HQ  or  is  true  using 
binary  antipodal  signaling  in  white,  Gaussian  noise.  The 
well  known  probability  of  error,  P(E),  for  this  case  is 

P(E)  =  erfc(V2E/NQ)  (11) 

where  the  message  bits  are  assumed  equally  likely  and  where 
the  error  function  compliment  (erfc)  is  defined  as 

i  00  o 

erfc(x)  =  -  S  exp(-y  /2)  dy  (12) 

72tt  x 

Thus,  in  the  case  of  interference  by  white,  Gaussian  noise, 
the  performance  of  the  SS  system  is  identical  to  a  conven¬ 
tional  system. 

The  real  advantage  of  the  use  of  an  SS  system  is  in 
the  presence  of  a  strong  narrowband  interfering  signal 
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centered  about  the  carrier  frequency.  Let  the  interfering 
signal,  n.(t),  be  zero  mean  and  Gaussian  with  PSD  ST(w)  as 
shown  in  Figure  3«  The  received  signal  is  now 


r(t)  =  s ( t )  +  n(t)  +  nj(t)  (13) 

and  hence 

r'(t)  =  -V2E/T  m(t)  cos(w0t)  +  c(t)  [n(t)  +  nj(t)] 

=  V2E/T  m(t)  cos((j)0t)  +  n"(t)  (14) 

The  autocorrelation  function  of  n"(t)  is  given  by 

R^Vtg)  =  ECn-C^)  n"(t2)} 

=  Etct^)  c(t2)} 

•E{[n(t1)  +  n^t-^]  [n(t2)  +  nj( tg )]} 

=  Rc^tl_t2^  CE£n( -bq )  n(t2)}  +  E{n(t1)  nj(t2)} 

+  E[nj(  t^)  n(t2)}  +  ECnjt^)  nj(t2)}]  (15) 

Since  n(t)  and  rij(t)  are  independent  of  each  other,  the  auto¬ 
correlation  becomes 

R"<Vt2)  =  Rc(tl-t2)  ^N(/2)  *(  VV  +  WV^  (16) 

With  T  =  t^-t2,  "^e  autocorrelation  of  the  interfering  signal 
is  found  by  taking  the  inverse  Fourier  transform  of  S  T( ci) ) . 

Hj(t)  =  ^  J*  Sj(co)  do) 

— OO 

*1  (i)n+B/2 

=  -J-  S  0  ,  (Jq/2B)  eja)T  dw 

W0-B/2  V 

(  •  Jo  —  (V)  (17) 
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Figure  3.  Power  Spectral  Density 
of  Interfering  Signal,  nj(t). 

The  autocorrelation  function  of  n"(t)  then  becomes  the 
followings 

Jn  sin(-rrBT) 

R"(T)  =  RC(T)  (Nq/2)  6  (T )  +  RC(T)  — -gT -  COS(C00T) 

Jn  sin(-rrBT) 

=  (Nq/2)  6(t)  +  Rc(t)  -  -  cos((jjQT )  (18) 


Now  let  B-*0.  This  corresponds  to  narrowband  jamming. 
Then  the  PSD  of  R"(t)  is  the  followings 

S"(<u)  =  J*  (Nq/2)  6 (t )  e“^WT  dT 

-oo 

oo 

+  J1  rc(t)  J0  cos(&)0t)  e“JWT  dT 

mmOO 

sin2[(co-d)0)Tc/2] 


=  (Nq/2)  +  (1/2)  J0T 


0  [(w-w0)t c/2y 


(19) 


+ 


(i/2)  JqTc 


sin2[(tO+W0)Tc/2] 

[(«+W0)Tc/2]2 


S"(w)  is  shown  in  Figure  4.  Since  the  data  pulse  duration, 

T,  is  typically  several  orders  of  magnitude  greater  than  the 
code  pulse  duration,  T  ,  the  bandwidth  of  c(t)*nT(t)  is  much 
wider  than  the  bandwidth  of  c(t)*s(t).  Over  the  message 
bandwidth,  the  PSD  of  c(t)*nj(t)  is  approximately  flat  and  is 
viewed  by  the  correlation  detector  as  an  additive  "white" 
noise  term.  Using  this  approximation,  equation  (19)  becomes 

S"(w)  «  (Nq/2)  +  (J0Tc/2)  (20) 

The  probability  of  error  in  deciding  whether  m(t)  was  a  plus 
one  or  a  minus  one  is  thus 

P(E)  -  erfc[V2E/(N0+J0Tc)]  (21) 

If  the  PN  code  was  not  used  in  the  system,  correspond¬ 
ing  to  c(t)  =  1,  then  the  autocorrelation  of  the  noise  is 
(with  B-*0) 

R"(r)  =  (Nq/2)  6(t)  +  JQ  cos(w0t)  (22) 

and  the  PSD  of  the  noise  is 

S"(w)  »  (Nq/2)  +  (Jq/2)  6(w+(D0)  +  (Jq/2)  6(u-W0)  (23) 

To  determine  the  associated  P(E)  for  this  case,  the  mean  and 
variance  of  the  sufficient  statistic,  e,  must  be  calculated. 
In  this  case 

T 

e  *  2E/T  X  m(t)  cos2((j)Qt)  dt 
0 
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Figure  4.  Power  Spectral  Density  of  n" (t). 


+  -J2E/1  S  n(t)  cos(cont)  dt 
0  u 

T 


+  V2E/T  S  nT(t)  cos(wnt)  dt 
0  J  o 


(24) 


The  conditional  means  are  given  by 


E{e|Hn}  =  2E/T  S  E{m(t)|Hn}  cos2(cont)  dt 


+  •JZE/T  S  E{n( t )  Ihq}  cos(w0t)  dt 


+  V2E/T  S  ECrij(t)|H0}  cos ( (i)q t )  dt 


=  2E/T  J*  cos^(a)nt)  dt  =  E  =  en 
0  u  u 


(25) 


and 


EtetH^  *  -E  *  e1 


(26) 
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The  conditional  variance  of  the  statistic,  given  either 
hypothesis,  is  identical  and  is  given  by  the  following 
(Ref.  20sl77)s 

VAR(e[HQ)  »  VARtelH^ 

T  T 

=  2E/T  X  X  R"(t,u)  cos(ti)nt)  cos(ti)nu)  dt  du 
0  0  u  u 

=  (ENq/2)  +  (EJ0T/2)  =  d2  (27) 

Thus,  the  probability  density  function  (PDF)  of  the  suffi¬ 
cient  statistic,  e,  given  either  HQ  or  H1  is  true,  is 
Gaussian  with  mean  eQ  or  e-^,  and  variance  d  .  The  associated 
probability  of  error  in  estimating  the  message  sent  is  then 

P(E)  =  P(e>0lH1)  P(H1)  +  P(e<0[H0)  P(HQ) 

=  P(e<0lH0) 

r  _  2  • 

1  1  ®0^  T  /  9TT 

=  -z=  x  exp  4  - dx  =  erfc  L  g  (28) 

V2 rrd  -«>  (_  d  cT  J  LVN0  J0. 

Equations  (21)  and  (28)  can  now  be  compared  to  see  the 
effect  of  the  code,  c(t).  If  NQ  is  considered  negligible, 
then  the  effect  of  the  narrowband  interference  on  the  system 
performance  was  reduced  from  JQT  to  Jqtc*  The  improvement 
ratio  is 

J0T  T 

jV  *  5  “  N  (29) 

J0Tc  Tc 

This  improvement  thru  the  use  of  the  PN  code  is  commonly  re¬ 
ferred  to  as  the  processing  gain  of  the  SS  system.  Notice 
that  it  is  roughly  equivalent  to  the  ratio  of  the  SS  trans- 
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mission  bandwidth  to  the  message  bandwidth.  Thus,  it  is  de¬ 
sirable  to  use  an  SS  system  in  an  environment  where  strong 
narrowband  interference  is  encountered.  It  should  be  noted 
that  this  analysis  assumes  that  the  jamming  signal  is  uncor¬ 
related  with  the  code,  c(t).  For  a  repeater  jammer,  the  re¬ 
sult  of  this  analysis  would  not  be  correct  because  of  the 
correlation  between  the  code  and  jamming  signal. 

The  fundamental  aspects  of  SS  communications  systems 
that  have  been  presented  in  this  chapter  should  equip  the 
reader  sufficiently  for  the  development  that  follows.  In  the 
next  chapter,  a  general  waveform  equation  will  be  presented. 
This  equation  will  be  used  as  the  primary  tool  in  the  analy¬ 
sis  of  the  proposed  intercept  receiver  model  of  chapter  IV. 

In  addition,  two  generic  SS  modulation  formats  will  be  intro¬ 
duced  and  shown  to  be  equivalent  to  the  general  waveform 
equation. 
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Ill  Generic  Modulation  Formats 

In  this  section,  a  waveform  equation  will  he  intro¬ 
duced  in  a  general  form  that  encompasses  many  different  mod¬ 
ulation  formats  that  are  already  well  described  in  the  liter¬ 
ature.  Two  generic  SS  modulation  formats,  format  A  and  B, 
and  their  associated  transmitters  will  be  introduced  that 
will  be  analyzed  in  a  later  section  for  their  vulnerability 
to  intercept.  Finally,  it  will  be  shown  that  the  two  generic 
modulation  formats  fit  the  general  waveform  equation. 

General  Waveform  Equation 

The  SS  modulation  formats  will  be  modeled  in  quadra¬ 
ture  component  form  as  (Ref.  2): 

s(t)  =  I(t)  cos[co0t  +  9(t)] 

+  Q(t)  sin[(oQt  +  6(t)]  (30) 

where 

I(t)  =  E  anp(t-nTc)  (3D 

Q(t)  =  E  V(t_nTc“Tc/2)  (32) 

and 

p(t)  =  pulse  shaping  function, 

a  =  code  and/or  data  dependent  term, 
n  constant  on  the  interval 
nTc-V2  <  t  <  nTc+Tc/2, 

b  *  code  and/or  data  dependent  term, 
n  constant  on  the  interval 
nTc  <  t  <  (n+l)Tc, 

I 
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Tq  =  code  pulse  duration, 

and 

6(t)  =  phase  term  caused  by  data 
modulation  (format  B  only). 

The  pulse  shaping  function,  p(t),  is  determined  by 
the  specific  RF  modulation  technique  used.  Three  methods  of 
RF  modulation  will  be  considered  in  this  thesis*,  offset  qua¬ 
ternary  phase  shift  keying  (O-QPSK);  minimum  shift  keying 
(MSK);  and  sinusoidal  frequency  shift  keying  (SFSK).  These 
methods  are  considered  for  the  following  reasons  (Ref.  5)* 

1)  These  techniques  exhibit  minimal  power  spectral 
occupancy,  which  is  desirable  to  an  SS  system  designer  who 
wants  to  minimize  detectability  of  the  signals. 

2)  Modern  system  non-linearities  and  power  efficien¬ 
cies  constrain  the  modulation  format  to  a  constant  envelope 
type,  which  these  are.  This  has  the  added  benefit  of  deny¬ 
ing  a  potential  interceptor  any  amplitude  modulation  infor¬ 
mation. 

For  these  reasons,  the  SS  system  designer  would  very  likely 
use  one  of  the  three  techniques!  O-QPSK,  MSK,  or  SFSK. 

For  O-QPSK  modulation,  the  pulse  shaping  function  is 
given  by 

1  -T  /2  <  t  <  T  /Z 

P(t)  -  -  ^  ^ 

0  elsewhere  (33) 


For  MSK,  the  pulse  shaping  function  is  given  by 


p(t)  = 


cos 


rrt 


-T q/2  <  t  <  Tc/2 

u 

0  elsewhere  (3^) 

For  SFSK,  the  pulse  shaping  function  is  (Ref.  1) 

-Tq/2  <  t  <  Tc/2 

elsewhere  (35) 


P(t)  = 


cos 


[nt  1  .  /Wt\ 

Tc ' 5  “TO 


The  three  pulse  shapes  are  shown  in  Figure  5.  These  three 
forms  of  modulation  are  identical  except  for  the  pulse  shap¬ 
ing  function  (Refs.  1,2,5,11,  and  18). 

To  get  a  clear  picture  of  the  general  waveform  de¬ 
scribed  by  equations  (30)  thru  (32),  and  to  understand  the 
relationships  between  the  terms,  consider  E  p(t-nTc)  and 

E  p(t-nT  -T  /2)  for  0-QPSK  pulse  shaping.  These  sums  are 
n  c  c 

shown  in  Figure  6.  The  second  pulse  stream  is  offset  from 
the  first  by  T  /2  (half  a  code  bit).  If  the  following  values 

C 

of  an  and  bn  are  assumed,  then  I(t)  and  Q(t)  are  as  shown  in 
Figure  6: 

a_2  =  1  a_^  =  -1  aQ  =  -1  a1  =  1  a2  =  _1 

b_2  =  1  b_i  =  !  bQ  =  -1  b2  =  -1  b2  =  1 


Notice  that  over  the  interval  0  <  t  <  T  .  there  are 

c 

two  terms  of  I(t)  involved  and  one  term  of  Q(t).  Thus,  to 
describe  the  waveform  over  any  interval  T  seconds  long  re- 
quires  two  equations.  The  representation  for  the  nth  inter¬ 
val,  CnT  <  t  <  (n+l)T  ],  is  the  followingt 
c  c 
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Figure  5»  Pulse  Shapes  for  O-QPSK,  MSK,  and  SFSK  Modulation 

s(t)  =  anp(t-nTc)  cos[wQt  +  ©(t)] 

+  bnp( t-nTc-Tc/2)  sin[u0t  +  ©(t)] 


for  nT  <  t  <  nT  + 
c  c 


V2 

=  an+1p[t-(n+l)Tc]  cos[ci)0t  +  ©(t)] 

+  bnp[t-nTc-T(J/2]  sin[w0t  +  ©(t)] 
for  nTc  +  1^2  <  t  <  (n+l)Tc 


(36) 
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I 

- — -sHL*. 


This  representation  is  general  enough  to  accomodate  a  wide 
variety  of  digitized  voice  modulation  techniques  and  will  be 
used  throughout  this  thesis,  adapted  to  fit  both  generic  mod¬ 
ulation  formats  considered.  Note  that  the  only  difference  in 
the  two  parts  of  equation  (36)  is  that  in  the  second  expres¬ 
sion,  the  term  an+1  has  replaced  an,  and  p[t-(n+l)Tcl  has 
replaced  p(t-riT  ).  Because  this  difference  involves  only  a 

C 

substitution  of  terms,  only  the  first  equation  will  be  used 
in  the  intercept  receiver  calculations  in  chapter  IV, 

The  two  generic  modulation  formats  and  their  trans¬ 
mitter  models  will  next  be  described.  The  primary  differ¬ 
ence  in  the  two  models  is  that  one  scheme  codes  the  data 
prior  to  modulation  on  carriers  in  quadrature  (format  A), 
while  the  other  scheme  directly  modulates  a  carrier  which  is 
then  coded  at  RF  on  quadratures  (format  B). 

Generic  Modulation  Format  A 

Figure  7  illustrates  the  block  diagram  for  the  trans¬ 
mitter  of  this  model.  Pulse  width  modulation  (PWM)  is  used 
to  develop  a  "digitized”  baseband  signal.  The  baseband 
signal  is  split  into  an  in-phase  and  quadrature  channel  where 
each  is  coded  by  one  of  two  orthogonal^  PN  sequences,  c^(t) 
and  Cg(t). 

In  this  model,  a  normalized  speech-voltage  waveform, 
m(t),  is  sampled  every  T  seconds  to  produce  a  stream  of  voice 
samples,  m^  =  m(kT).  The  kth  sample  is  pulse  width  modulated 

^Two  codes  are  orthogonal  if  the  expected  value  of 
their  product  is  zero,  i.e.,  E(cj( t)c^( t)}  =  0. 
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Figure  7.  Generic  Transmitter  Model  A 


to  generate  a  pulse  of  variable  width  on  the  interval 
kT  <  t  <  (k+l)T,  as  is  illustrated  in  Figure  8.  Typical 
examples  are  shown  for  positive,  negative,  and  zero  level 
samples.  The  pulse  width,  rfc,  for  the  kth  voice  sample  is 
(Ref.  4) 


Tk  =  (1+mk)T/2  (37) 

Each  T  second  sample  interval  of  the  pulse  width  mod¬ 
ulated  data  stream,  d(t),  is  divided  into  two  T/2  sections. 

The  first  section  is  routed  to  the  in-phase  channel  and  the 
second  to  the  quadrature  channel.  Each  section,  denoted 
dj(t)  and  d^(t),  contains  a  pulse  whose  width  is  doubled  and 
the  DC  level  is  shifted  to  zero  as  shown  in  Figure  8.  The 
data  streams  dj(t)  and  ^(t)  are  bipolar  with  value 

dI(t)  =  ±1  dQ(t)  =  ±1  df ( t)  =  dq( t)  =  1  (38) 

Note  that  there  is  a  natural  T/2  second  offset  in  the 

start  of  each  T  second  interval  of  dj(t)  and  d^(t).  For 
mathematical  convenience,  this  offset  will  be  ignored  real¬ 
izing  that  it  must  be  accounted  for  in  a  real  receiver.  Thus, 
it  will  be  assumed  that  each  interval  of  d^(t)  and  d^t)  is 
coincident  prior  to  multiplication  by  a  PN  code. 

The  I  and  Q  channel  data  are  modulated  by  two  orthog¬ 
onal  PN  codes,  Cj(t)  and  c^t),  with  value 

c-^t)  =*  ±1  Cq( t )  =  ±1  Cj(t)  =  c*(t)  =  1  (39) 

The  transitions  in  the  data  streams,  dj(t)  and  dqtt),  are 
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Figure  8.  Transmitter  Waveforms,  (a)  Sampled  voice,  m 
(b)  Pulse  width  modulated  samples,  d(t).  (c)  In-phase 

channel  data,  (d)  Quadrature  channel  data. 


quantized  such  that  they  are  synchronized  with  the  code 

transitions  so  that  the  data  is  completely  hidden  by  the  code. 

The  modulated  data  streams  are  next  multiplied  by  pulse 

shaping  functions,  I  p(t-nT  )  and  £  p( t-nT  -T  /2) ,  where  p(t) 

n  c  n  c  c 

is  as  defined  by  one  of  equations  (33)  thru  (35).  The  resul¬ 
tant  products  modulate  in-phase  and  quadrature  carriers,  and 
are  added  to  form  the  transmitted  signal 

s(t)  =  Z  cIndInP('t-nTc)  cos(Ci)Qt) 

+  l  cQndQnp(t-nTc-Tc/2)  sin(V}  (40) 

where 

"th 

cIn  =  n  term  of  Cj(t),  constant  on  the 

interval  nT  -T  /2  <  t  <  nT  +  T  /2, 
c  c  c  c 

CQn  =  nth  term  of  c^(t),  constant  on  the 
interval  nT  <  t  <  (n+l)T  , 

djn  =  value  of  dj(t),  constant  on  the 

interval  nT  -  T  /Z  <  t  <  nT  +  T  /2, 
c  c  c  c 

d^n  =  value  of  dq(t),  constant  on  the 

interval  nT  <  t  <  (n+l)T  . 

c  c 

For  the  nth  interval,  nT  <  t  <  (n+l)T  ,  the  output  signal  is 

c  c 

s(t)  =  cIndInp^ t_nTc^  c°8(w0t) 

+  cQndQnp(t"nTc_T(/2)  sin(wot) 
for  nTQ  <  t  <  nTc  +  TQ/2 

=  CI  ( n+1  )dI  ( n+1  )p^“ ( n+1  )Tc'1  cos(w0t) 

+  cQndQnp(t‘nTc‘Tc/2)  8in(«ot) 

for  nTc  +  1^/2  <  t  <  (n+l)Tc  (41) 
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Thus,  it  is  seen  that  the  output  of  transmitter  model  A  is 
identical  to  the  general  equation,  equation  (36),  with 

an  =  cIndIn’  bn  =  cQndQn’  and  an+l  =  cI(n+l)dI(n+l) * 

Generic  Modulation  Format  B 

In  the  transmitter  model  for  this  mqdulation  format, 
shown  in  Figure  9»  a  speech  voltage  waveform,  m(t),  is 
sampled  every  T  seconds.  The  kth  sample,  m(kT)  =  m^,  is 
mapped  to  a  phase  value,  ©k,  by  a  non-linear  differential 
amplitude  to  phase  rule.  The  relationship  between  the  voice 
sample  and  its  corresponding  phase  value  is 


e 


k 


f(mk)  +  e 


k-1 


(42) 


where  f )  is  the  non-linear  mapping  rule  and  ©k  1  is  the 
phase  value  for  the  previous  sample,  ®k_1.  The  mapping  rule 
is 


f(mk)  =  2  Cos 


(43) 


Thus,  mk  is  found  from  successive  phase  values  by 


(44) 


The  phase  values  directly  modulate  a  carrier,  cos((o0t),  to 
produce,  over  all  time,  cos[(i)Qt  +  0(t)3.  A  quadrature  car¬ 
rier,  sin[(i)Qt  +  0(t)],  is  derived  from  this  by  passing  thru  a 
90°  phase  shifter.  These  carriers  are  modulated  by  orthogonal 
PN  codes;  Cj(t)  and  c^(t),  and  pulse  shaping  data  streams, 

E  p(t-nTQ)  and  E  p( t-nTc-Tc/2) ,  where  these  terms  are  as 
described  previously.  The  in-phase  and  quadrature  waveforms 
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Figure  9.  Generic  Transmitter  Model  B 


are  added  to  give  the  transmitted  waveform 

s(t)  =  £  CjnPC't-nTc)  cos[w0t  +  0(t)3 

+  n  cQnP^t-nTc"Tc/2^  slnCw0t  +  (**5) 

For  the  interval,  nT  <  t  <  (n+l)T  ,  the  transmitted 

c  c 

signal  is 

s(t)  =  cInP('t-nTc)  cos[to0t  +  0n] 

+  cQnp(t-nTc‘Tc/2)  sinCw0t  +  ®n3 

for  nT  <  t  <  nT  +  T  /2 
c  c  c 

=  cI(n+l)p^t-(n+1)Tc^  cosCw0t  +  en] 

+  cQnpl-t'nTc-Tc//2^  sinC^o*  +  6n-^ 

for  nT  +  T  /2  <  t  <  (n+l)T  (46) 

c  c  c 

This  is  identical  to  the  general  equation  (36)  with  cIn  =  an 

and  cQn  ’  V 

It  has  been  shown  that  both  generic  formats,  A  and  B, 
fit  the  general  waveform  model  given  by  equation  (36).  A 
proposed  intercept  receiver  will  be  introduced  in  the  next 
section,  and  the  general  waveform  equation  will  be  used  to 
determine  its  outputs.  The  results  will  be  applied  to  gen¬ 
eric  formats  A  and  B. 


In  this  chapter,  a  proposed  intercept  receiver  model 


will  be  introduced.  Using  the  general  waveform  equation 
developed  in  Chapter  III  as  the  transmitted  waveform,  the 
receiver  outputs  will  be  calculated  and  examined  for  exploi¬ 
tation  and  potential  information  recovery.  It  will  be  shown 
that  under  appropriate  conditions,  this  receiver  can  success¬ 
fully  decode  the  transmitted  message  for  both  generic  wave¬ 
form  models  previously  introduced. 

Intercept  Receiver  Model 

The  proposed  intercept  receiver  model  is  shown  in 
Figure  10.  The  primary  components  are  the  sections  labeled 
the  Carrier  Reconstruction  Loop,  the  Format  A  Decoder,  and 
the  Format  B  Decoder.  Each  section  will  be  treated  in  de¬ 
tail  in  the  receiver  analysis  in  the  next  subsection. 

Some  motivation  for  this  particular  structure  is  in 
order.  This  configuration  represents  the  block  diagram  im¬ 
plementation  of  mathematical  massaging  of  the  general  wave¬ 
form  equation  to  isolate  modulation  terms  from  PN  code  terms 
without  the  aid  of  PN  code  correlation.  Recall  from  the 
discussion  of  the  properties  of  PN  codes  that  multiplying 
the  code  times  a  coherent  replica  eliminates  the  code,  that 
is,  c  (t)  a  1.  The  identical  effect  is  produced  by  proces¬ 
sing  thru  a  non-linear  element.  The  quadrupler  in  the 
Carrier  Reconstruction  Loop  and  the  square  law  elements  in 
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the  Format  A  Decoder  accomplish  this  purpose.  The  remain¬ 
ing  hardware  serves  only  to  isolate  the  uncoded  modulation 
terms  and  reduce  interference  and  noise.  The  purpose  of 
each  element  will  become  obvious  as  the  receiver  is  analyzed. 

Receiver  Analysis 

The  received  waveform  is  given  by 

r(t)  =  A(t)  Z  anp(t-nTc)  cos[w0t  +  cp(t)] 

+  A(t)  E  b  p(  t-nT  -T  /2 )  sin[wnt  +  cp(t)] 

n  n  C  C  v 

+  n(t)  (4?) 

where  A(t)  is  the  amplitude  of  the  received  signal;  cp(t)  is 
an  instantaneous  phase  term  that  includes  arrival  phase  off¬ 
sets,  Doppler  phase  shifts,  and  may  include  modulation 
(format  B);  n(t)  is  modeled  as  zero-mean,  white,  Gaussian 
noise  with  two  sided  PSD  height  Nq/2;  and  a^,  bR,  p(t-nTc), 
and  (t-nTc-Tc/2)  are  as  defined  in  equations  (31)  and  (32). 

Over  one  bit  interval  of  the  PN  code,  nT  <  t  <  (n+l)T  , 

c  c 

the  received  signal  is 

r(t)  =  Aanp(t-nTQ)  cos[w0t  +  cp(t)] 

+  Abnp(  t-nTc-Tc/2)  sin[(i)0t  +  cp(  t ) D  +  n(t) 

for  nTc  <  t  <  nTc  +  1^2  (48) 

and 

r(t)  =  Aan+1p[t-(n+l)Tc]  cos[ci)Qt  +  <p(t)] 

+  Abnp( t-nT(,-Tc/2)  sin[(i)0t  +  <p(t)]  +  n(t) 

for  nTc  +  1^2  <  t  <  (n+l)Tc  (49) 

In  equations  (48)  and  (49),  it  has  been  assumed  that  the 
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signal  amplitude  remains  constant  over  the  code  bit  inter¬ 
val.  Notice  that  the  only  difference  in  the  two  equations 
is  that  in  the  second,  aR+1  replaces  an  and  p[t-(n+l)Tc] 
replaces  p(t-nT  ).  Because  this  difference  involves  only  a 

V 

substitution  of  terms,  equation  (48)  will  be  used  as  the 

basis  for  the  succeeding  development. 

The  received  signal  is  filtered  by  a  bandpass  filter 

(BPF)  with  bandwidth  (BW)  B  proportional  to  the  reciprocal 

of  the  code  pulse  duration  (B  «  2/T  ).  Thus,  the  assumption 

is  being  made  that  the  transmission  bandwidth  B  includes  only 
4 

the  main  lobe  of  the  spectrum  of  the  transmitted  signal. 

The  filter  is  centered  at  the  carrier  frequency,  fQ.  The 
signal  out  of  the  filter  is  the  following: 

x(t)  =  [Aanp(t-nTc)  +  nQ(t)3  cos[a)0t  +  cp(t)] 

+  [Abnp(  t-nTc-Tc/2)  -  ng(t)3  sin[wQt  +  cp(t)]  (50) 

In  this  equation,  n  (t)  and  n  (t)  are  the  in-phase  and  quad- 

c  s 

rature  components  of  the  bandlimited  noise  at  the  output  of 
the  filter  (Ref.  22:237)*  The  noise  components  are  Gaussian 
with 

E{nc(t)}  =  E{ng(t)}  =  E{n(t)}  =  0  (51) 


E[n^(t)}  =  E[ng( t )}  =■  E[n  (t)}  =  (N0)*(Noise  BW)  (52) 

where  the  noise  bandwidth  equals  the  bandwidth  of  the  filter 
preceeding  it. 


^Refer  to  Figure  2.c  of  Chapter  II  for  a  typical 
baseband  spectrum  after  coding.  The  "main  lobe"  refers  to 
the  portion  of  the  spectrum  between  -2rr/T„  and  2rr/T  . 

W  C 
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At  this  point  the  filtered  signal,  x(t),  is  proces¬ 
sed  thru  the  Carrier  Reconstruction  Loop  and  the  Format  A 
Decoder.  The  Carrier  Reconstruction  Loop  will  be  considered 
first,  followed  by  the  Format  A  Decoder.  The  Format  B  De¬ 
coder  will  be  considered  last. 

Carrier  Reconstruction  Loop.  The  Carrier  Reconstruc¬ 
tion  Loop,  Figure  11,  serves  two  purposes  in  the  intercept 
receiver.  As  its  name  implies,  one  purpose  is  to  establish 
a  coherent  carrier  reference  in  the  receiver.  The  configur- 
ation  shown  is  a  special  case  of  the  N  power  loop  describ¬ 
ed  by  Lindsey  and  Simon  (Ref.  13:71-?4).  They  show  that  a 
coherent  carrier  is  produced  at  the  output  of  the  phase  lock 
loop  (PLL)  for  rectangular  shaped  data.  Since  the  operation 
of  the  PLL  is  well  documented  in  the  literature  (Ref.  8,13, 
19  and  21),  no  treatment  will  be  given  here.  It  will  be 
shown  in  Appendix  A,  however,  that  for  MSK  and  SFSK  pulse 
shaping,  a  signal  similar  to  the  rectangular  pulse  shape 
case  is  available  to  the  PLL  for  tracking.  Thus,  the  Carri¬ 
er  Reconstruction  Loop  will  generate  a  coherent  carrier  in 
all  cases  of  interest  here. 

The  second  purpose  that  the  Carrier  Reconstruction 
Loop  serves  is  to  provide  a  signal  to  the  Format  B  Decoder 
that  contains  the  desired  modulation  terms  but  no  code  dis¬ 
torted  terms.  This  will  become  apparent  in  the  analysis  of 
the  Format  B  Decoder. 

An  important  element  of  the  Carrier  Reconstruction 
Loop  is  the  quadrupling  circuit  described  in  Appendix  A. 
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Figure  11.  Carrier  Reconstruction  Loop. 

It  is  shown  in  the  appendix  that  the  quadrupling  circuit 
produces  a  signal  at  its  output  given  by 

■4  .  4  1 

v(t)  =  p4(t-nTc)  +  |  p\t-nTc-Tc/2) 

-  ^A4p2(t-nTc)  p2(t-nTc-Tc/2)j 
•  cos[4(jj0t  +  4cp(  t  )1 

[4 

I  anV3(t_nTc)  P(t-nTc-Tc/2) 

‘  I  anbnp(t_nTc)  P3(t-nTc-Tc/2)J 
•  sin[4(j}0t  +  4cp(t)] 

+  Nc(t)  cos[4o)0t  +  4cp(t)] 

+  Ns(t)  sin[4o)0t  +  4cp(t)] 


(53) 


where  N  (t)  and  N  (t)  are  noise  terms  defined  in  Appendix  A. 
c  s 

If  the  "signal"  component  of  v(t)  is  defined  to  be  the  terms 
near  not  multiplied  by  quadrature  noise  terms,  then  the 
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"signal"  component  at  the  output  of  the  quadrupling  circuit 
for  the  cases  of  O-QPSK,  MSK,  and  SFSK  pulse  shaping  are 
given  by 

4 

S  O-QPSK  =  •  I  cos[4co0t  +  Mt)3  (5^) 

SMSK  cos[4d30t  +  4p(t)]  (55) 

SSFSK  [^4  +  ^a4j  ( 1  )J  cos[4a)0t  +  4p(t)]  (56) 

where  J^(l)  is  the  Bessel  function  of  the  first  kind  of  order 
1  and  with  argument  1.  Thus,  regardless  of  the  pulse  shaping 
used  by  the  transmitter,  for  the  generic  modulation  formats 
considered  here  the  output  of  the  quadrupling  circuit  can 
always  be  expressed  as 

v(t)  =  -V  cos[4(i)0t  +  4cp(t)] 

+  Nqv t)  cos[4w0t  +  4cp(t)] 

+  Ns(t)  sin[4d)0t  +  4p(t)]  (57) 

This  result  is  divided  in  frequency  by  four  to  give 

v*(t)  =  [-V  +  Nc(t)]  cos[w0t  +  tp(t)] 

+  Ng(t)  sin[w0t  +  cp( t>D  (58) 

This  signal,  then,  is  presented  to  the  Format  B  Decoder  and 
the  PLL.  As  is  demonstrated  in  Appendix  A,  the  signal  to 
noise  ratio  (SNR)  is  typically  quite  large  here,  and  hence 
the  noise  contributions  are  negligible.  Then  for  all  cases 
of  interest,  a  trackable  signal  is  presented  to  the  PLL  and 
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coherent  reference  signals  at  the  carrier  frequency  can  he 
produced.  Also,  the  Format  B  Decoder  is  given  a  signal 
containing  the  modulation  terms  (in  the  case  of  modulation 
format  B)  that  is  free  of  PN  code  terms  and  significant 
noise  terms. 

With  coherent  carrier  reference  signals  available, 
the  intercept  ‘receiver  is  able  to  decode  transmissions  from 
both  modulations  format  A  and  B.  The  first  of  these  to  be 
considered  will  be  format  A. 

Format  A  Decoder.  The  Format  A  Decoder  is  shown 
again  in  Figure  12.  Here,  the  signals  in  the  in-phase  and 
quadrature  channels  will  be  calculated.  Outputs  will  be 
calculated  for  0-QPSK,  MSK,  and  SFSK  pulse  shaping  of  the 
data.  It  will  be  shown  that  when  only  one  transmission  is 
received,  message  decoding  is  not  possible  at  the  channel 
outputs.  However,  when  two  messages  from  different  trans¬ 
mitters  are  received  simultaneously,  it  will  be  shown  that  it 
is  possible  to  decode  both  messages  if  the  transmissions  ar¬ 
rive  at  the  receiver  within  the  same  PN  code  bit  interval 
(i.e.,  the  transmitters  are  synchronized  and  approximately 
equal  distances  from  the  receiver). 

When  one  transmission  is  received  the  input  to  the  de¬ 
coder  is  given  by  equation  (50),  repeated  here. 

x(t)  «  [Aanp(t-nTQ)  +  nQ(t)]  cos[(i)0t  +  Cp( t )3 

+  [Abnp(  t-nTc-Tc/2)  -  ng(t)]  sin[w0t  +  qp(t)]  (50) 
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Figure  12.  Format  A  Decoder 


This  signal  is  mixed  with  coherent  reference  signals  in  the 
in-phase  and  quadrature  channels  to  produce 

xI(t)  =  lanp(t-nTc)  +  lnc(t) 

+  [fanp^t-nTc^  +  Inc(t)]  cosC2W0t  +  2cp(t)] 

+  [lbnP^t"nTc“Tc//2^  "  sinC20)0t  +  2cp(t)3(59) 

and 

xQ(t)  =  |hnp(t-nTc-Tc/2)  -  |ns(t) 

-  [^nP( ‘t-nTc-Tc/ 2 )  -  |ng(t)]  cos[2d)0t  +  2cp(t)] 

+  [fanp^t-nTc^  +  K(t)]  sinC2^0t  +  2cp ( t ) D  (60) 

These  signals  are  filtered  thru  ideal  low  pass  filters  (LPF) 
with  BW  =  B,  the  transmission  bandwidth.  Because  the 
carrier  frequency  is  much  larger  than  the  transmission  band¬ 
width,  only  the  first  two  terms  in  each  equation  are  passed. 
The  output  of  the  filters  is  then: 

x*(t)  =  |anp(t-nT0)  +  |ne(t)  (61) 

x*(t)  »  -  K(t)  <62) 

The  signals  are  now  processed  thru  square  law  devices  to 
.  > 

produce  signal  terras  independent  of  coding.  Recalling  that 
2  2 

=  b^  =  1,  the  output  of  the  square  law  devices  is: 

2 

Uj(t)  =  |  p2(t-nTc)  +  £n2(t) 

+  |anP(t-nTc)nc(t)  (63) 
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(64) 


2 

“Qtt)  ■  5  P2(t-nTe-T,/2)  » 

-  fbnp(t-nTc-Tc/2)ns(t) 

At  this  point,  it  is  necessary  to  include  the  specific 
pulse  shaping  functions.  Using  equations  (33),  (34)  and  (35) 
in  (63)  and  (64),  the  output  of  the  square  law  device  for 
O-QPSK,  MSK,  and  SFSK  pulse  shaping  is  as  follows.  For 
O-QPSK : 


V*5  =  5  +  5no(t)  +  fv'c(t) 

(65) 

Vt5  =  K(t>  '  Kns(t) 

(66) 

For  MSK  pulse  shapings 

„  /  .X  _  A2.  A2  /2Trt\ 

ul(t)  "  H  F  \%) 

+  |an  cos^^nc(t)  +  Jn2(t)  (67) 

u  f  +  %  _  A2  A2  /2rrt\ 
uQ(t)  ■  h  ■  s 

-  |bn  sin^ns(t)  +  £n2(t)  (68) 

The  output  of  the  quadrupler  for  SFSK  pulse  shaping  is  more 
difficult  to  evaluate,  so  intermediate  steps  will  be  includ¬ 
ed  here.  In  this  case 

Uj(t)  =  £2cos2plE  .  1 
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+  fan  cos[?;  -  i  sin(^)]nc(t)  +  3nc(t) 

_  A2  .  A2  P 2trt  1  .  /W\l 

-  8  +  8  C0S[^;  -  2  slnVT;;J 

,  A  fret  1  .  /VrrtVl  ,  .*  .  1  2, 

+  2an  cos[f;  -  5  sin\wJ  c(  +  5nc(t) 

=  f 2  +  Jn(l/2)  cos^(l-2n)^j 

+  fan  g  Jn(l/4)  cos[(l-4n)^jnc(t)  +  Jn*(t)  (69) 

xy. 

where  Jn(|i)  is  the  n  order  Bessel  function  of  the  first 
kind  with  argument  p,  and  where  (Ref.  22s  116) 

v  cos[wQt  +  (d)sin(ojmt]  =  v  E  Jn(u)  cos[(d)c+  nwm)t]  (?0) 

The  output  of  the  square  law  device  in  the  quadrature  channel 
for  SFSK  pulse  shaping  is,  similarly 

uq( t )  =  §2-  |2S  Jn(l/2)  cos[\l-2n)^j 

-  fbn  £  Jn(l/4)  8in|(l-4n)^Jn8(t)  +  Jn*(t)  (71) 

Before  stating  the  outputs  of  the  Format  A  Decoder, 
it  is  necessary  to  consider  the  effect  of  the  final  LPF. 

This  ideal  filter  has  a  bandwidth  W  on  the  order  of  the  re¬ 
ciprocal  of  the  sampling  rate  at  the  transmitter  (W  «  l/T), 
the  information  bandwidth.  If  the  PSD  and  autocorrelation 
function  of  the  output  were  of  interest,  then  it  would  be 
appropriate  to  compute  these  characteristics  at  the  input  to 
this  filter  at  this  time  and  derive  the  equivalent  expres¬ 
sions  for  the  output.  However,  determining  the  output  in 


^5 


( 


the  time  domain  from  the  output  autocorrelation  function  is 
not  possible.  For  example,  two  PN  codes  have  the  same  auto¬ 
correlation  function,  but  their  time  history  can  be  complete¬ 
ly  different.  Since  a  time  domain  expression  for  the  output 
is  of  primary  importance  here,  the  following  representation 
be  used  to  describe  the  portion  of  the  signal  passed  by  the 
ideal  low  pass  filter:  [  ]Lp. 

The  output  of  the  Format  A  Decoder  for  the  different 
pulse  shaping  functions  when  one  transmission  is  received, 
then,  is  as  follows.  For  O-QPSK: 


For  MSK: 


yI('t)  = 
= 


yI(t)  = 

yQ(t)  = 


5nc(t)  +  fannc(t) 
ys{t)  -  fbnns(t) 


5nc(t> 


LP 

LP 


LP 

(72) 

LP 

(73) 

(74) 

(75) 


For  SFSK: 


(76) 

(77) 

Recall  now  that  for  modulation  format  A,  a„  and  b, 

n  n 

are  composed  of  modulation  terms  and  code  terms.  Notice  in 
the  output,  equations  (72)  thru  (77)»  that  there  are  only 
noise  terms  or  coded  modulation  terms  multiplied  by  noise. 
Thus,  it  is  clear  that  for  the  case  where  only  one  trans¬ 
mission  is  received,  no  information  can  be  recovered  at  the 


yj(t)  = 

yQ(t)  = 


*nc(t) 

'blit) 


LP 

LP 
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output.  However,  it  will  now  "be  shown  that  for  the  case 

where  two  simultaneous  transmissions  are  received  by  the 

interceptor  that  arrive  within  the  same  PN  code  bit  interval, 

message  reconstruction  is  possible. 

Consider  the  received  waveform  for  two  simultaneous 

transmissions.  A  PN  code  bit  offset  of  At  <  TV  will  be  as- 

c 

sumed  in  the  analysis.  A  carrier  phase  difference  of  will 

be  assumed  between  transmitter  1  and  transmitter  2,  also. 

The  received  waveform  for  the  interval  nT  <  t  <  nT„  +  T  /Z 

c  c  c 

is 

r(t)  =  A1anlp(t-nT(J)  cos((i)Qt) 

+  AibniP(t-nTc-Tc/2)  sin((i)Qt) 

+  A2an2P( t-nT  -At)  cos(w0t  +  $) 

+  A2bn2p( t-nT  -Tc/2-At)  sin(wQt  +  0) 

+  n(t)  (78) 

where  the  subscripts  1  and  2  refer  to  the  contributions  from 
transmitters  1  and  2.  With  a  =  cos  $  and  3  =  sin  $,  the 
signal  at  the  output  of  the  front  end  filter  is 

x(t)  =  [A1anlp(t-nTc)  +  oA2an2p^ t~nTc  "  ^ 

+  3A2bn2p( t-nTc-Tc/2-At)]  cos(w0t) 

*  C*!  bnip(t-nTc-Tc/2)  +  aA2bn2p(t-nTc-Tc/2-At) 

-  3A2an2p( t-nTc-At)  -  ng(t)]  sin(ojQt)  (79) 

The  calculations  necessary  to  compute  the  outputs  become 
quite  lengthy  and  complicated  at  this  point,  hence,  they 
have  been  included  in  Appendix  B.  The  outputs  of  the  Format 
A  Decoder  for  two  simultaneous  transmissions  are  listed  as 
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follows.  For  0-QPSK  pulse  shaping* 


A2anlan2 


yi(t)  2AlA2anlan2  +  [|A1 

?]lp  +  [|n 


5Li.A2a  b 

2  2n2  n2 


:(t) 


LP 

LP 


a„,  +  0A^a._„  +  ^A„a__„^  nc(t)J 


llanl  T  2A2an2  +  lA2an2 


LP 


yQ(t)  2AlA2bnlbn2 


[lAlA2an2bnl]LP 

-  +  [5ns(t)]LP 

-  (Kbnl  +  I 


"A2bn2  +  2A2an2 


:)nS(t)]LP 


For  MSK  pulse  shaping,  the  outputs  are: 


yI(t)  =  §AlA2anlan2 

-[!AlA2anlbn2  Sin(^)]Lp  +  [Ve^LP 
yQ(t)  "  SAlA2bnlbn2  C°S(^t) 

+[lAlA2an2bnl  Sin(^)]Lp  +  [K™]lP 


(80) 


(81) 


(82) 


(83) 


The  corresponding  terms  for  SFSK  pulse  shaping  ares 
yI(t)  =  |A1A2anlan2  cos^^J2(l/4) 

‘  [lAlA2an2bnl  ain(^)  J0(1/,4)jLp  +  [inc(t)]LP  (84) 
yQ ^  =  5AlA2bnlbn2  C°S(^^)  Jq^A) 

+  [!AlVn2bnl  8in(%;)  J0(l/4)]Lp  +  frs^LP  (85) 
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Before  proceeding  further  in  the  developement,  it  is 
instructive  to  examine  the  output  equations  in  detail.  First 
consider  what  the  interceptor  is  looking  for:  terms  that 
contain  the  information  but  are  free  of  PN  code  terms  and  are 
not  multiplied  by  noise.  For  modulation  format  A,  recall 
that  =  cIndInl.  an2  =  cIndjn2,  V  -  cQndQnl,  and 

^n2  ~  cQndQn2  w^ere  the  "c"  terms  are  PN  code  terms,  the 
"d"  terms  are  the  desired  data  or  information,  and  where  the 
numeric  subscripts  refer  to  transmitter  1  and  2  [see  defin¬ 
itions  under  equation  (40)3.  These  terms  appear  in  the  out¬ 
put  equations  in  the  following  combinations: 


anlan2 

dInldIn2 

(86) 

anlbn2 

cIncQndInldQn2 

(87) 

an2bnl 

cIncQndIn2dQnl 

(88) 

an2bn2 

cIncQndIn2dQn2 

(89) 

bnlbn2  = 

dQnldQn2 

(90) 

The  code-free  terms  are  thus  the  ones  that  contain  either 
anl*an2  or  b^ • bn2 .  Then  the  first  term  of  each  of  the 
Format  A  Decoder  output  equations,  equations  (80)  thru  (85)» 
is  the  desired  term  from  which  the  interceptor  would  try  to 
reconstruct  the  message. 

Also  notice  in  the  output  equations  that  the  trans¬ 
mitter's  relative  carrier  phase  offsets  and  the  PN  code  bit 
offset  play  a  major  role  in  the  relative  amplitudes  of  the 


49 


terms.  The  terms  containing  code  have  a  PSD  that  is  low 
compared  to  the  terms  without  code,  so  it  can  be  argued 
heuristically  that  they  remain  at  a  relatively  low  and 
approximately  constant  power.  Thus,  the  amplitude  of  the 
first  term  is  of  primary  importance  in  its  detection. 

Consider  now  the  ideal  case  where  the  carriers  arrive 
in  phase,  the  PN  code  bits  exactly  coincide,  and  the  noise 
power  is  negligible.  This  corresponds  to  a  =  1,  0=0, 

At  =  0,  and  n  (t)  =  n  (t)  =  0.  The  decoder  outputs  under 
these  ideal  conditions  are  given  in  Table  I.  These  same  re¬ 
sults  are  duplicated  in  Table  II  except  that  the  substitu¬ 
tions  for  the  "a"  and  "b"  terms  have  been  made,  and  the  "n" 
subscript  has  been  dropped  since  it  only  indicates  that  the 
results  are  for  the  nth  interval.  Examination  of  Table  II 
shows  that  regardless  of  the  pulse  shaping  used,  the  output 
of  the  receiver  contains  the  essential  speech  ingredients 
multiplied  by  a  constant. 

How,  then,  is  the  message  reconstructed  from  these 
outputs?  First,  it  is  assumed  that  the  outputs  are  normal¬ 
ized  by  some  gain  control  circuity  so  that  the  output  in 
each  case  is 

yjtt)  =  dndI2  (91) 

and 

yQ(t)  =  dQ1dQ2  (92) 

Keep  in  mind  that  filtering  the  output  had  the  effect  of 
integrating  over  the  sample  period,  T.  The  above  equation 

I 
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TABLE  I 


Format  A  Decoder  Output 
Under  Ideal  Conditions. 
General  Waveform  Case. 


Pulse 

Shaping 

yx(t) 

yQ(t) 

O-QPSK 

2AlA2anlan2 

tAlA2bnlbn2 

MSK 

5AlA2anlan2 

5AlA2bnlbn2 

SFSK 

5AlA2bnlbn2JS<1A> 

TABLE  II 

Format  A  Decoder  Output 

Under  Ideal  Conditions. 

Modulation  Format  A  Case. 

Pulse 

Shaping 

yj(t) 

yQ(t) 

O-QPSK 

iAlA2dIldI2 

2AlA2dQldQ2 

MSK 

5AlA2dIldI2 

5AlA2dQldQ2 

SFSK 

5AlA2dQldQ2  J0  ^ ^ 

is  the  output  over  ove  PN  code  chip  interval,  T  .  The  out- 

c 

put  over  all  time  is 

yI(t)  =  di:L(t)dI2(t)  (93) 

yQ(t)=  dQl(t)dQ2(t)  (94) 

Let  d^(t)  and  d2(t)  he  the  waveforms  at  the  output  of  the 
pulse  width  modulators  of  transmitters  1  and  2  respectively 
(see  Figure  7).  Since  d1(t)  and  d2(t)  can  over  any  sample 
interval  represent  a  positive,  negative,  or  zero  sample, 


then  there  are 

nine  cases 

which  must  be 

considered. 

1) 

d1(t) 

d2(t) 

positive 

positive 

2) 

d-^t) 

d2(t> 

positive 

negative 

3) 

d1(t) 

d2(t) 

negative 

positive 

4) 

d^(t) 

d2(t) 

negative 

negative 

5) 

dx(t) 

d2(t) 

positive 

zero 

6) 

dx(t) 

d2(t) 

negative 

zero 

7) 

dx(t) 

d2(t) 

zero 

positive 

8) 

dx(t) 

d2(t) 

zero 

negative 

9) 

dx(t) 

d2(t) 

zero 

zero 

Each  of  the  nine  cases  is  illustrated  in  Figure  13. a 
thru  13. i.  Consider  case  2,  Figure  13.b,  where  d^(t)  repre¬ 
sents  a  positive  sample  and  d2(t)  represents  a  negative 
sample.^  Here,  and  in  all  nine  cases,  and  r2  represent 

^The  reader  is  encouraged  to  review  the  Generic 
Modulation  Format  A  section  of  chapter  III  for  clarifica¬ 
tion  of  this  development. 
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Figure  13. a.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  1. 
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T  o-T/2 


Figure  13.c.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  3. 
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Figure  13 .d.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  4. . 
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Figure  13. e.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  5* 
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Figure  13.f.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  6. 
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Figure  13. h.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  8. 


Figure  13* i.  Data  Signals  and 
Format  A  Decoder  Outputs  - 
Case  9. 
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the  pulse  widths  of  d^(t)  and  d2(t)  respectively  over  one 
sample  interval  of  duration  T.  Notice  that  d-^t)  has  been 
divided  into  two  sections  of  duration  T/2.  The  first  sec¬ 
tion  has  a  constant  value  over  the  interval,  so  it  can  be 
represented  by  a  pulse  of  width  T/2.  This  pulse  width  is 
doubled  and  routed  to  the  in-phase  channel  at  the  transmit¬ 
ter  to  produce  the  signal  d-^t).  Similarly,  the  second 
section  of  d^(t)  can  be  represented  by  a  pulse  of  width 
r^-T/2.  This  also  is  doubled  in  width  and  routed  to  the 
quadrature  phase  channel  of  the  transmitter  to  produce  the 
signal  d^^t).  The  signals  dI2(t)  and  dQ2(t)  are  derived 
from  d2(t)  in  identical  fashion.  The  figure  also  shows  the 
outputs  of  the  receiver,  yj(t)  and  y^t),  as  given  by  equa¬ 
tions  (93)  and  (94).  By  a  similar  treatment,  the  outputs 
are  developed  for  the  remaining  eight  cases. 

From  equation  (37) »  it  is  clear  that  the  pulse  widths, 
and  t2,  of  the  message  sample  values,  mkl  and  mk2,  are 
related  by  the  following: 

T1  =  (1+mkl)T/2  (95) 

T2  =  (l+mk2)T/2  (96) 

Then  the  message  sample  values  are: 

mkl  =  2rl/T  ‘  1  (97) 

ate  =  2t g/T  -  1  (98) 

Thus,  to  recover  the  message  samples,  and  r2  must  be  de¬ 
termined  from  the  outputs  of  the  receiver,  y^(t)  and  yg(t). 
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One  way  that  the  interceptor  might  proceed  is  to  try 
to  recover  m^  and  independently  of  each  other,  so  he 
must  come  up  with  some  decoding  rule  to  accomplish  this. 
Notice  that  for  any  particular  case,  the  output  may  contain 
zero,  one,  or  two  zero  crossings  in  the  sample  interval. 

The  times  to  the  first  and  second  zero  crossing  (if  any)  in 
yj(t)  are  designated  as  t-[1  and  tI2,  respectively.  In  yQ(t), 
these  are  designated  as  t^  and  tQ2.  ^  should  be  clear  that 
the  interceptor  can  not  know  whether  a  particular  zero  cross¬ 
ing  is  caused  by  transmitter  1  or  2.  Therefore  he  must  make 
some  assumption.  He  assumes,  rightly  or  wrongly,  that  when 
only  one  crossing  occurs  in  y^t),  it  is  caused  by  transmit¬ 
ter  1,  and  when  only  one  occurs  in  y^(t)  that  it  is  caused 
by  transmitter  2.  When  two  zero  crossings  occur  in  y^(t), 
the  first  is  assumed  to  be  caused  by  transmitter  1  and  the 
second  by  transmitter  2.  In  yQ ( i: ) ,  the  first  is  assumed 
caused  by  transmitter  2  and  the  second  by  transmitter  1. 

To  summarize: 

tj1  -*  transmitter  1  ( ) 
tjg  "*■  transmitter  2  (f2) 
tQ!  ■*  transmitter  2  (r2) 

*Q2  "*  ‘transmi't'ter  1  (T^) 

Under  the  above  assumptions,  the  following  estimates  of  T, 

/  -L 

and  r2  can  be  made.  From  the  in-phase  channel  (yj(t))i 

=  tj^/2  (99) 

T2  =  tI2/2  (100) 
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(101) 


From  the  quadrature  channel  (yQ(t)): 

^1  =  ^Q2+  TV2 

T2  =  (tQ1+  T)/2  (102) 

Using  these  estimates  in  equations  (97)  and  (98)  gives  the 
following  decoding  rules.  For  the  in-phase  channel  (y-j-(t)): 

~  -  1  (103) 

**2  =  tI2/T  -  1  (104) 

For  the  quadrature  channel  (yq(t))i 

“kl  =  tQ2/T  U°5> 

mk2  =  tQ2/T  (106) 

When  only  one  transition  occurs  in  either  channel,  the  rule 
is 

^kl  =  tIl/T  -  1  (!°7) 

for  the  in-phase  channel  and 

mk2  =  tQ1/T  -  1  (108) 

for  the  quadrature  channel.  The  obvious  flaw  with  this 
method  is  that  on  the  average,  this  will  result  in  a  correct 
decision  being  made  only  fifty  percent  of  the  time  since  the 
interceptor  can't  know  which  transmitter  caused  which  tran¬ 
sition. 

Since  the  human  ear  does  a  good  job  of  discriminating 
simultaneous  voices,  a  better  scheme  would  be  to  determine 
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In  this  case 


the  sum  of  the  message  samples,  m^  +  m^. 

“kl  +  mk2  =  2( fi+Tg )/T  -  2 


(109) 


The  problem  now  becomes  one  of  estimating  instead  of 
and  r2  individually. 

Consider  Case  1,  Figure  13. a,  where  d^(t)  and  d2(t) 
both  represent  positive  samples.  At  yQ(t),  suppose  that  a 
wrong  assumption  is  made:  that  tQ^  is  caused  by  transmit¬ 
ter  2.  Mathematically,  the  assumption  is: 


Then 


*1  =  (tQ2+  T)/2 


T2  =  (tQl+  T)/2 


^1  +  ^2  ^tQl+tQ2^//2  +  T 


(110) 


(111) 


(112) 


In  reality,  t^  was  caused  by  transmitter  1.  Mathematically, 


this  is: 


Then 


~  ("tQi+  T)/2 

t2  ~  (V+  ^^2 

T1  +  T2  =  (tQl+tQ2)/2  +  T  =  TX  +  f; 


(113) 


(11*0 


(115) 


Thus,  even  though  the  wrong  assumption  was  made,  by  adding 
the  result,  the  sum  of  the  samples  can  be  reconstructed  un¬ 
ambiguously  on  each  sampling  interval.  This  result  occurs 
for  each  of  the  nine  possible  cases.  In  each  case,  the  sum 
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of  the  two  samples  can  be  recovered  without  having  to  make 
an  assumption  about  which  transmitter  caused  a  particular 
transition  at  the  output. 

When  the  preceeding  analysis  is  performed  for  the 
case  of  three  simultaneous  transmissions,  the  results  are 
not  as  successful.  In  this  case,  the  object  is  to  detect 
three  zero  crossings  at  the  Format  A  Decoder  output  per  sam¬ 
ple  period,  T.  However,  in  all  cases,  only  two  of  the  three 
expected  zero  crossings  occur.  Since  the  interceptor  does 
not  know  which  two  transmitters  cause  the  zero  crossings 
that  he  does  detect,  he  must  guess  according  to  some  prede¬ 
termined  rule.  Regardless  of  the  rule,  one  third  of  the  in¬ 
formation  has  been  lost  and  performance  is  degraded  over  the 
two  transmission  case. 

Thus,  for  modulation  format  A,  intercept  and  decoding 
of  the  message  for  the  case  of  a  single  transmission  is  not 
possible  without  knowledge  of  the  PN  code  used.  However,  as 
has  been  shown,  for  two  simultaneous  transmissions  and  under 
ideal  conditions,  the  modulation  format  A  waveform  can  be 
intercepted  and  decoded  without  knowledge  of  the  PN  codes 
used.  As  was  pointed  out,  the  fidelity  of  the  reconstructed 
message  will  be  degraded  by  additive  noise,  the  phase  offset 
between  the  transmitter  carriers,  and  the  offset  between  the 
code  bits  at  the  receiver.  For  three  simultaneous  trans¬ 
missions,  only  partial  decoding  is  possible  since  one  third 
of  the  information  is  lost.  In  the  next  section,  the  Format 
B  Decoder  will  be  analyzed.  It  will  be  shown  that  decoding 
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a  modulation  format  B  waveform  does  not  require  more  than 
one  transmission  as  is  true  for  the  format  A  case. 


Format  B  Decoder.  Transmissions  from  generic  trans¬ 
mitter  model  B,  Figure  9,  are  processed  by  the  Format  B  De¬ 
coder,  shown  in  Figure  14.  This  configuration  is  a  logical 
choice  for  the  interceptor  if  he  assumes  the  following: 

1)  The  phase  information  is  differentially  encoded. 

2)  The  rule  used  to  map  the  voice  samples  to  phase 
values  is  linear. 

The  second  assumption  is  not  ludicrous  when  one  considers 
that  for  any  continuous  mapping  rule,  small  variations 
around  some  operating  point  are  nearly  linear.  In  any  case, 
the  interceptor  must  assume  some  mapping  rule,  and  a  linear 
rule  is  as  attractive  as  any.  Based  upon  these  assumptions, 
the  interceptors  estimate  of  the  encoding  rule  is 


6k  =  (Klmk  +  K2>  +  6k-l 


(116) 


where  and  K2  are  constants,  m^  is  the  value  of  the  sample 
on  the  k  sample  interval,  and  and  are  the  encoded 

phase  values  for  the  kth  and  (k-l)st  samples.  Thus,  the 
interceptors  decoding  rule  is  the  following: 


m 


k 


K? 

ek-l)  ‘ 


(117) 


The  interceptors  first  task,  then,  is  to  detect  the  phase 
values,  and  6^^.  Then,  he  must  use  his  decoding  rule, 
equation  (117).  These  tasks  are  implemented  by  the  hard¬ 
ware  of  Figure  14. 
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Assuming  a  negligible  noise  level,  the  input  signal 
v  (t)  is  multiplied  by  a  sinusoid  at  the  carrier  frequency 
to  give 

f(t)  =  -V  cos[w0t  +  cp( t ) H  sin[ti)Qt  +  cp0] 

=  sin[2w0t  +  cp(t)  +  cp0]  +  |  sinQp(t)  -  cpQ]  (118) 

where  cpQ  is  a  phase  offset  between  the  sinusoid  and  the  in¬ 
put  signal,  and  where  cp(t)  is  a  phase  term  composed  of  encod¬ 
ed  phase  terms,  © ( t ) ,  and  some  Doppler  phase,  ©d(t),  associ¬ 
ated  with  the  relative  motion  between  the  transmitter  and 
intercept  receiver.  Thus 

cp(t)  =  ©(  t)  +  ed(t)  (119) 

Equation  (118)  is  filtered  by  a  lowpass  filter  to  give  the 
following: 

g(t)  =  sin[cp(t)  -  cp0]  (120) 

This  is  processed  to  give 

h(t)  =  Sin-1{sin[cp(t)  -  cpQ]}  =  Cp( t )  -  cpQ  (121) 

If  it  is  assumed  that  the  Doppler  phase  varies  slowly 
enough  so  that  it  can  be  considered  constant  over  a  sampling 
interval  and  approximately  equal  from  one  interval  to  the 

S  "t 

next,  then  the  following  results  occur.  For  the  (k-1) 
interval 

Mt)  =  hk_1  =  ek_1  +  ed  -  <pQ  (122) 
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For  the  interval 

h(t)  =  hk  =  ek  +  ed  .  cpQ  (123) 

Then 

P(t)  =  hk'-  hk-1  =  ek  -  «k_x  (124) 

The  interceptor  has  now  reconstructed  the  differential 
phase  and  will  process  it  according  to  his  decoding  rule, 
equation  (117).  The  output  is  then 


%  -  -  W  -  § 


(125) 


Recall  that  the  transmitter  used  the  following  non¬ 
linear  encoding  rule: 

\  - 2  Cos'1(^?)  +  \-i  <126> 

The  function 

-  »)E_1  =  2  Cos_1(i^)  (127) 

is  graphed  in  Figure  15  over  the  range  -1  ^  mk  s  1.  Using 
equation  (12?)  in  (125),  the  sample  estimate  at  the  output 
becomes 


A 


(128) 


As  a  measure  of  the  accuracy  of  this  estimate,  it  is 
useful  to  compute  the  Mean  Square  Error  (MSE).  This  is 
defined  as 
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Figure  15.  Format  B  Encoding  Rule,  PDF  of  Message  Samples 
and  Interceptors  Estimated  Encoding  Rule. 


MSE  =  E{(mk-mk)2} 

=  ECra2}  -  2E(mkmk}  +  E{m2} 

'  ECm^}  -  ^  E[mk  Cos-1^)] 

+  E«  *  ^2  E  [cos-1^)] 

'  ^  E[Cos'1(Tk)]  +  E  j|]  <12?» 

To  evaluate  the  MSE  further  requires  knowledge  of  the  proba 
bility  density  function  of  mk.  Recall  that  the  parameter 
mk  is  an  amplitude  sample  of  a  speech  voltage  waveform,  m(t). 
An  adequate  model  for  a  speech  voltage  waveform  distribution 
is  given  by  the  Laplace  distribution  (Ref.  19:46).  The 
associated  probability  density  function  is  given  by  the  fol¬ 
lowing  (Ref.  14): 

p(m)  =  %  exp(-y|m|)  Y  >  0  (130) 

The  parameter  y  is  determined  as  follows.  Since  the  function 
Cos  ^[(l+mk)/23  is  not  defined  for  |mkl>l,  the  original  voice 
waveform  was  preprocessed  such  that  the  probability  of  m. 

K 

exceeding  1  would  be  very  small.  If  mk  is  arbitrarily  con- 
■ trained  such  that 

Pr[|mk!>l]  £  0.01  (131) 

•»  T  ilt  be  such  that 

X1  \  exp(-Ylml)  dm  >  0.99  (132) 

-1 
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Solving  this  for  v  gives 


Y  >  -  ln( 0,01)  k  4.6052  (I33) 

Using  this  value,  the  probability  density  function  of  m^, 
also  plotted  in  Figure  15»  is  given  by  the  following*. 

p(m)  =  2.3026  exp(-4.6052tml)  (134) 

Its  first  and  second  moments  are 


E(mk}  =  0 

(135) 

eU^I  =  -2  =  0.0943 

(136) 

Using  this  probability  density  function,  the  MSE, 
equation  (129),  can  now  be  evaluated.  The  expectations  in¬ 
volving  Cos-1[(l+mk)/2]  were  computed  numerically  to  give 
the  following  equation  for  the  MSE: 

MSE  =  0.0943  +  (0.20)/^ 

+  (4.53  -  4.20K2  +  K2)/kJ  (137) 

The  minimum  MSE  is  found  by  taking  the  partial  derivatives 
of  the  MSE  with  respect  to  and  K2,  set  them  equal  to  zero, 
and  solve  simultaneously  for  and  K2*  This  procedure  gives 
K-^  =  -1.2  and  K2  =  2.1.  The  minimum  MSE  is  now  computed  to 
be 


MSEmin  =  °-011 


(138) 


The  minimum  MSE  normalized  to  the  a  priori  voice  variance  is 

MSE  • 

- =  0.II63  (139) 

ECnC} 
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The  interceptors  linear  decoding  rule,  equation  (127), 
has  also  been  plotted  in  Figure  15  using  the  values  of  K1  and 
K2  computed  above.  Notice  how  closely  this  approximates  the 
actual  non-linear  rule  used  for  small  variations  about  zero. 
Thus,  under  the  ideal  conditions  assumed,  the  output  of  the 
Format  B  Decoder  closely  approximates  the  transmitted  mes¬ 
sage. 

Unfortunately,  ideal  conditions  are  seldom  encount¬ 
ered.  The  Doppler  phase  can  cause  significant  errors.  As 
a  worst  case  example  for  a  fixed  intercept  receiver,  con¬ 
sider  the  case  of  the  transmitter  aboard  an  aircraft  moving 
at  Mach  2  relative  to  the  receiver.  Assume  a  transmitter 
frequency  of  fQ  =  250  MHz.  The  Doppler  frequency,  fd,  is 
(Ref.  13) 

fd  =  f0(l-v/c)  (140) 

where  v  is  the  aircraft  velocity  (662.72  m/sec)  and  c  is  the 

O 

speed  of  light  (2.998  x  10  m/sec).  Thus,  for  the  worst  case 

fd  =  249.9994474  MHz  (141) 

For  a  signal  of  the  form  cos[(i)Qt  +  ©d(t)],  the  instan¬ 

taneous  phase  is  Wgt  +  ©d(t).  The  instantaneous  angular 
frequency  is 

to  =  |g[w0t  +  ©d(t)]  (142) 

Thus,  the  Doppler  angular  frequency,  cod  isi 

0)d  =  2nfd  =  <o0  +  |t[ed(t)] 
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(143) 


“d -  2"f0 * 

Then,  the  Doppler  phase  is 

t 

6d(t)  =  2t r-J*  (frf-fQ)  dt 

u  0  u 

=  -2rr(552.6)t  (144) 

Suppose  that  the  sampling  rate  is  f  =  l/T  =  8  kHz.  With 
t  =  T,  the  Doppler  phase  shift  over  one  sample  interval  is 

ed(t)  =  -2tr(552.6)T 

=  -0.43  radians/interval 
=  -24.87  degrees/interval  (145) 

This  is  a  significant  error.  Under  most  cases,  the  Doppler 
phase  must  be  handled.  If  it  is  assumed  that  ©d(t)  changes 
at  approximately  a  constant  rate  (corresponding  to  a  constant 
aircraft  velocity),  then  the  following  results  occur.  For 
the  (k-l)st  intervals 

hk_i  =  ®k_i  +  ed[(k-l)T  +  t]  -  cp0  (146) 

For  the  k^*1  intervals 

hk  =  \  +  ed(kT  +  t)  -  <pQ  (147) 

Thus 

pk  =  Vhk-1  =  (V®k-l}  +  ted(kT+t)  -  edC(k-l)T+t]}  (148) 
and 

K 

z(t)  =  g  (®k“®k-l>  "  %  +  i1ted(kT+t)  ■  ®dC(k-l)T+t]}  (149) 
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It  is  desired  that  the  last  term  in  equation  (149)  be 
zero  to  eliminate  Doppler  phase  errors.  One  way  to  accom¬ 
plish  this  is  to  estimate  the  Doppler  phase  and  subtract  it 
thru  some  type  of  decision  directed  feedback  circuitry. 

This  is  also  referred  to  as  data  aided  feedback,  decision 
feedback,  modulation  removal,  or  inverse  modulation  (Ref. 
8:219-225  and  13:64).  The  interested  reader  is  referred  to 

these  sources  for  a  comprehensive  look  at  this  technique. 

« 

In  this  section,  it  was  shown  that  the  Format  B  De¬ 
coder  can  decode  the  intercepted  signal  if  the  additive 
noise  is  considered  negligible,  if  the  Doppler  phase  can  be 
estimated,  and  if  the  rule  used  at  the  transmitter  to  map 
amplitude  samples  of  the  message  into  phase  values  is  nearly 
linear,  at  least  for  small  variations  about  zero.  Unlike 
the  case  of  modulation  format  A,  decoding  modulation  format  B 
does  not  require  two  simultaneous  transmissions. 


mammam 
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V  Conclusions  and  Recommendations 

Conclusions 

This  thesis  has  shown  that  for  a  class  of  modulation 
formats,  intercept  and  message  decoding  of  pseudo-noise  en¬ 
coded  spread  spectrum  waveforms  is  possible  without  know-  • 
ledge  or  duplication  of  the  spreading  code  used,  under  appro¬ 
priate  conditions.  The  essential  element  necessary  to  accom¬ 
plish  this  was  processing  of  the  received  signals  through 
non-linear  elements  (square  law  devices  and  fourth  law  de¬ 
vices).  The  non-linear  devices  had  the  effect  of  producing 
terms  that  had  been  stripped  of  code  terms. 

A  general  waveform  equation  was  introduced  as  the 
primary  tool  in  the  analysis.  This  waveform  equation  could 
be  adapted  to  fit  a  number  of  polyphase  shift  keyed  modula¬ 
tion  schemes.  It  was  demonstrated  that  this  general  wave¬ 
form  equation  encompassed  O-QPSK,  MSK,  and  SFSK  modulation 
schemes  when  the  appropriate  pulse  shaping  function  was  used. 

Two  generic  modulation  formats,  designated  as  format  A 
and  format  B,  and  their  associated  transmitter  models  were 
introduced.  It  was  shown  that  these  formats  fit  the  general 
waveform  equation.  In  format  A,  the  message  was  split  into 
quadratures  and  encoded  at  baseband  by  two  orthogonal  PN 
codes  prior  to  RF  modulation.  In  format  B,  quadrature  car¬ 
riers  were  directly  phase  modulated  with  identical  informa¬ 
tion  prior  to  encoding  by  two  orthogonal  PN  codes  at  RF. 
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An  intercept  receiver  model  was  proposed.  The  general 
waveform  equation  was  used  as  the  input  to  the  receiver  and 
the  receiver  signals  were  computed.  The  receiver  was  divided 
into  three  main  sections:  the  Carrier  Reconstruction  Loop; 
the  Format  A  Decoder;  and  the  Format  B  Decoder. 

The  Carrier  Reconstruction  Loop,  composed  of  a  quad¬ 
rupling  circuit  and  a  PLL,  was  analyzed  first.  It  was  shown 
that  regardless  of  the  format  used  (A  or  B)  or  the  pulse 
shaping  used  on  the  data  (0-QPSK,  MSK,  or  SFSK),  there  was 
a  term  produced  at  the  output  of  the  quadrupler  that  did  not 
contain  code  terms.  This  term  could  be  tracked  by  the  PLL 
to  produce  coherent  reference  signals  at  the  carrier  frequen¬ 
cy.  In  addition,  for  format  B,  this  term  contained  the  de¬ 
sired  message  information. 

It  was  shown  that  for  modulation  format  A,  decoding 
the  message  without  knowledge  of  the  PN  code  was  not  possible 
for  this  receiver  when  only  one  transmission  was  received. 

It  was  demonstrated  that  when  two  simultaneous  transmissions 
were  received,  message  reconstruction  of  both  transmissions 
was  possible  under  the  appropriate  conditions.  The  output  of 
the  Format  A  Decoder  was  calculated  for  this  case  and  the 
results  showed  that  successful  decoding  depends  upon  the  dif¬ 
ferences  in  the  two  transmitters  carrier  phases,  the  arrival, 
offset  between  the  PN  code  bits  at  the  receiver,  and  additive 
noise.  When  ideal  conditions  were  assumed  (no  carrier  phase 
difference,  no  code  offset,  negligible  noise)  it  was  demon¬ 
strated  that  both  transmitted  messages  could  be  reconstructed. 
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Under  non-ideal  conditions,  the  performance  was  degraded. 

For  three  simultaneous  transmissions,  one  third  of  the  de¬ 
sired  information  was  not  present  at  the  output. 

For  modulation  format  B,  decoding  of  the  message  was 
possible  even  when  only  one  transmission  was  received.  Fac¬ 
tors  which  affected  the  success  of  decoding  this  format  were 
additive  noise,  Doppler  phase  shift  associated  with  the  rela¬ 
tive  motion  between  the  transmitter  and  receiver,  and  the 
extent  to  which  the  baseband  modulation  was  linear. 

Recommendations 

The  following  recommendations  for  further  study  are 
made  by  the  author.  First,  the  signal  to  noise  ratio  of  the 
quadrupling  circuit  output  should  be  calculated  for  the  case 
of  MSK  and  SFSK  modulation  to  verify  the  negligible  noise 
power  assumption  made  in  Appendix  A.  Secondly,  since  the 
results  in  this  thesis  were  developed  for  primarily  ideal 
conditions,  the  analysis  should  be  extended  to  quantify  the 
conditions  under  which  successful  intercept  and  decoding  is 
possible.  It  is  also  recommended  that  the  signals  and  receiver 
performance  be  simulated  on  a  computer  to  test  the  validity  of 
the  results  of  the  extended  analysis.  If  possible,  hardware 
tests  should  be  performed.  Finally,  it  is  recommended  that 
other  intercept  methods  be  investigated,  such  as  correlating 
a  received  signal  with  a  delayed  version  of  itself. 
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Appendix  A:  Quadrupling  Circuit  Calculations 

In  this  appendix,  the  output  of  the  quadrupling  cir¬ 
cuit  shown  in  Figure  16  will  he  calculated  for  input  signals 
of  the  form  of  equation  (48),  repeated  here. 

r(t)  =  Aanp(t-nTc)  cos[wQt  +  ep(t)] 

+  Abnp( t-nTc-Tc/2)  sin[co0t  +  cp(t)] 

+  n(t)  (48) 

The  output  equation  developed  here  will  he  shown  to  he  con¬ 
sistent  with  results  in  the  open  literature.  The  amplitudes 
of  the  signal  components  of  the  output,  v(t),  will  he  calcu¬ 
lated  for  the  cases  of  0-QPSK,  MSK,  and  SFSK  pulse  shaping. 
Finally,  the  signal  to  noise  ratio  at  the  output  will  he  cal¬ 
culated  for  O-QPSK  pulse  shaping. 

The  received  signal  is  filtered  to  give 

x(t)  =  [Aanpa  +  nc(t)]  cos[w0t  +  cp(t)] 

+  [Ahnpb  -  ns(t)]  sin[wQt  +  ep(t)]  (150) 

In  equation  (150),  p_  and  p,  have  been  substituted  for  nota- 
tional  compactness  for  p(t-nTc)  and  p( t-nTc~Tc/2)  respective¬ 
ly.  Also,  n(t)  has  been  represented  by  its  narrowband  quad¬ 
rature  components,  n  (t)  and  n  (t)  (Ref.  22»237)«  Again,  for 

c  s 

notational  compactness,  the  arguments  of  n  (t)  and  n_(t)  will 

c  s 

be  dropped.  Raising  equation  (150)  to  the  fourth  power  and 
keeping  only  those  terms  in  the  vicinity  of  4to0  gives  the 
output  of  this  circuit,  v(t). 
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v(t)  = 


Figure  16.  Quadrupling  Circuit. 

[§A>a  +  lA4pb  ‘  i*A4papb]  J 

+  [lA\iVapb  "  iA\ibnpapb]  ainCHt  + 

+  [^A3anp2  *  IA3anpap?)nc  +  (lA3Vlpb  ' 

+  'AVatb’Vs  +  (~2Aanpa)ncns 

*  (H*kS  *  -  fr24K 

*  QA2Pb  -  ^2paK  -  i  ncns  +  (I'VaK 

+  (‘lAVb)n3  +  5  nc  +  B  ns]  cosC4(i)0t  +  '4^<t)] 

+  [(lA3Vapb  -  ^VbX  +  (lA3anpapb  '  XvJK 

+  (lA2p?  -  lA2p!K"S  +  (_2Abnpb)ncns 
+  ("lAanpa)ncns  +  Q#ZaMfbK 

*  (*i*2a„Vapb)ns  +  (lAbnpb)"3  +  (lAanpa)ns 

*  Xns  '  !ncns]  sinC4'*J0t  +  M*)]  (1 


U51) 


83 


It  would  be  nice  to  compare  this  equation  with  a  known 
result  for  verification.  To  do  this,  consider  the  case  of 
O-QPSK  (rectangular)  pulse  shaping  and  no  PN  coding.  In  the 
equation,  this  corresponds  to  letting  an  =  t>n  =  1,  and 
Pa  =  p^  =  1.  With  rp( t )  =  0,  the  output  of  the  quadrupling 
circuit  is  then: 


v(t)  =  cos(4<o0t) 


* 

+  -A^n  +  A-^n  +  3A2n  n  -  ^An  n* 

C  S  C  S  £  C  £ 

I 

+  2Ancns  -  Kns  +  lAnc  -  iAns 

*  sno  +  sns]  cos(Ht) 

+  [a3"o  +  A\  -  2Ancns  -  lAncns 

+  ^  *  |An3  *  Vs 

+  5ncns  -  Kns]  8in(lK0ot) 


(152) 


An  identical  result  is  obtained  by  Butman  and  Lesh 


(Ref.  3)*  The  input  signal  in  their  article  is 

r(  t) '  =  s(  t ) '  +  n(  t) ' 

where  n(t)'  is  zero  mean,  white,  Gaussian  noise  and 


(153) 


s(t)'  =  J2A  cos(wnt  +  2itr/4) 


i  =  0,1, 2, 3  (154) 


Notice  that  s(t)’  can  take  on  any  one  of  four  values  desig¬ 
nated  Sq,  s£,  s,,,  and  s y  Figure  (17)  is  the  signal  space 
representation  of  r(t)'  in  the  r|,  r£  plane.  In  the  figure, 
n(t)'  is  shown  in  its  quadrature  component  representation, 
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Figure  17 .  Signal  Space  Representation 
of  the  Received  Signal,  r(t)’. 

n  (t)'  and  n  (t)'.  The  filtered  signal  is  given  as 
c  s 

r*(t)'  =  R(t)’  cos[a)0t  +  2irr/4  +  $>(t)] 


where 


1 

R( t ) '  =  CCV2A  +  nc(t)’]2  +  n2(t)'}2 


and 

$(t)'  =  Tan_1Cns(t)V[V2A  +  nc(t)']} 

xu 

Raising  this  to  the  4  power  and  retaining  only  terms 
around  4o)q  gives  the  output, 

v(t)'  =  R-^^'  cos[4<i)0t  +  44>(t)] 


(1 55) 

(156) 

(157) 

(158) 
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Using  equations  (1 56),  (157).  several  trigonometric  identi¬ 
ties,  and  dropping  the  arguments  of  nc(t)'  and  ns(t)'f  equa¬ 
tion  (158)  becomes  the  following: 

v(t)'  =  cos(4(i)Qt) 

+  [^AV  -  |A2n2'  ♦  |A2n2'  .  |V5An c'n2. 

*  H'4'  *  I-3-  ♦  ♦  $«*•]  oos(4»0t) 

*  [-VZA3ns.  +  3A2nc.ns"  -  ^An2'^' 

+  1An2'  -  K-2'  +  tn2'ns']  sin(4M0t)  (159) 

Equations  (152)  and  (159)  do  not  appear  to  be  identical  until 

the  differences  in  the  signal  representation  are  taken  into 

f 

«  account.  For  the  example  at  hand,  equation  (48)  can  be 

written  equivalently  as 

r(t)  =  s ( t )  +  n(t) 

=  V2 A  cos[(i)0t  +  (2i+l)rr/4] 

+  n(t)  i  =  0,1, 2, 3  (160) 

Figure  18  shows  the  corresponding  signal  space  representation 
for  this  case.  From  Figure  17  and  18,  it  is  clear  that 

nc(t)’  =  ^§nc(t)  -  ^§ns(t)  (161) 

and 

ns(t)’  =  ^|nc(t)  +  ^n8(t)  (162) 

When  the  differences  in  the  signal  and  noise  representations 
(  are  taken  into  account,  equation  (159)  can  be  shown  to  be 

equivalent  to  (152).  Thus,  equation  (151)  appears  correct. 
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Figure  18.  Signal  Space  Representation 
of  the  Received  Signal,  r(t). 

The  output,  equation  (151),  is  now  rewritten  as 
v(t)  =  [§A4p4  +  |a4p4  -  ^A4pjp^l  cos[4«0t  +  lKp(t)] 

+[tA\Vapb  *  §A\Vapb]  sinC^t  +  4,p<t>] 

+  Nc(t)  cos[4<i30t  +  4cp(t)] 

+  Ns(t)  sin[4<jJ0t  +  ^cp(t)]  (163) 

where  N  (t)  and  N  (t)  are  noise  terms  whose  definitions  are 
c  s 

apparent  from  a  comparison  of  equations  (151)  and  (I63). 

The  "signal"  terms  will  he  defined  to  he  the  first  two 
expressions,  and  the  noise  terms  the  last  two  expressions 


8? 


of  equation  (163).  The  signal  terms  will  now  be  evaluated 
for  the  cases  of  0-QPSK,  MSK,  and  SFSK  pulse  shaping. 

The  pulse  shaping  function  for  O-QPSK  is  given  by- 
equation  (33)  •  From  this,  it  is  clear  that 


pm(t-nTc)  = 


nT  -  T  /2  <  t  <  nT  +  T  /Z 
c  c  c  c 


elsewhere 


(164) 


pm ( t-nT  -T  / 2 )  = 

C  G 


1  nT  <  t  <  (n+l)T 
c  c 


elsewhere 


(165) 


These  expressions  are  substituted  into  equation  (163) 
signal  component  becomes  the  following: 


O-QPSK 


=  cos[4a)0t  +  4p(t)] 


(166) 


The  average  power  in  this  signal  is: 


T  4  2 

Ps  =  ^  <f  "I  cosC(i)0t  +  4{p(t)]  dt 

A®  8 

=  |  =  0.1250  A° 


(16?) 


From  equation  (34),  the  following  is  true  for  MSK 


modulation: 


pm(t-nTc)  =  (-1)"“  cosm(rrt/Tc) 


for  nTc  -  Tc/2  <  t  <  nTc  +  Tc/2 


elsewhere  (168) 
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pm(t-nT  -T  /2)  =  (-l)nmsinra(Trt/T  ) 

C  C 

for  nTc  <  t  <  (n+l)Tc 

=  0  elsewhere  (169) 

These  are  substituted  into  the  signal  terms  of  equation  (I63) 
to  gives 


3A  cos' 


•  cos[4  Qt  +  4  (t)] 


A* 

B  a 


„b 
n  n 


•  sin[4a)Qt  +  4cp(t)] 


cosC4oj0t 


+  4<p(t)3 


+ 


+ 


[Hb  "  I§Banbn]  C0S[Ht  + 

its  +  fl5anbn]  +  Mt)  - 


(170) 


Since  the  first  term  in  this  equation  is  the  only  one  near 
4(i)q •  it  is  the  only  one  that  actually  appears  in  the  output 
of  the  quadrupling  circuit.  The  last  two  terms  would  be  fil¬ 
tered.  Thus,  the  signal  output  for  MSK  modulation  is 


SMSK  =  ^5  cos[4o)0t  +  4<p(t)] 


(171) 


I 


f 

4 


and  its  average  power  is 


< 


PS  =  5?52  =  °'0011  A® 


(172) 


A  similar  treatment  can  be  given  for  SFSK  modulation. 
Using  equation  (35),  "the  following  equations  result: 

pm(  t-nT0 )  =  (-1)™  co5"fe  -  l  sin  *2*1 

*•  c  xc  J 

for  (nTc  -  Tc/2)  <  t  <  (nTQ  +  TQ/2) 

=  0  elsewhere  (173) 

pm(t-nTc-Tc/2)  =  (-I)™1  sinm^  -  £  sin 


for  nT  <  t  <  (n+l)T 
c  c 


=  o 


elsewhere  (17*0 


Substituting  these  expressions  into  equation  (163)  and  let- 
ting  q  =  ^  sin^p^  the  signal  term  for  SFSK  modulation  is: 

S  -  4/rrt  .  4  Art  \ 

ssfsk  ■  |b  cos  ^  "  7  s  sln  (t ~  -  y 

-  ^  cos2(^  -  q^  sin2^r  -  q^j  cos[4(o0t  +  4cp(t)] 

+  [l  Vn  oos3(t  '  l)  sin(^  -  l) 

*  5  anbn  oos (i;Vn3(^r  ■  4)j  3inCi|"ot  *  075) 

Using  trigonometric  substitutions,  rearranging  terms,  and 
substituting  for  q  gives 
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SSFSK  =  it  cosCZ,'wot  + 


[its  "  H5anbn]  cos[^0t  +  ^  -  sin 

nbn]  cosL^t  +  4cp(t)  -  ^  +  sin  ^1 

u  c  cj 


4  U 
+  £A_  .  ZPi 

1§S  l2Ba 


(176) 


Letting  oj  =  2tt/T  =  2trf  ,  this  becomes 

C  L  U 


4 

S  =  it  cosCi^ot  +  4cp(t)] 

+  [l28  "  128anbn]  cos^W0t  +  4cp(t)  +  2wct  -  sin(2wct)] 

sin[4o)0t  +  4cp(t)  -  2wct  +  sin(2wct)] 


+  +  2A^a  b  ] 

128  128  nbnj 


To  evaluate  this  further,  let 


ws  =  H  +  2wc 
wd  =  "  2wc 
%  =  2wc 


(177) 

(178) 

(179) 

(180) 


Then 


’SFSK  ~  64  cos^i<’Wo't  +  **$(*)  3 

+  [lI5  "  l§5anbn]  C0SEti)8*  +  V11)  “  sinCo^t)] 

I1  4  4  1 

128  *■  ffaVnH^v  +  4|*,(t)  +  sin'V,;1 


(181) 
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Since 


cos[wt  +  v  +  (i  sin(wmt)]  =  £  Jn(|i)  cos[(w+nwm)t  +  v](l82) 
then  equation  (181)  becomes  (with  all  substitutions  removed)! 


SSFSK  =  it  cos^4wot  +  M*)] 


+ 


+ 


i§5  '  I§8antn]  l  Jn(1)  008  jW  -(n-l)^£  +  ^t)] 

12B  +  lf8anbn]  n  Jn^  cos^4w0t  +  (n-l)^p  +  4cp(t)] 

C  (183) 


From  a  table  of  Bessel  functions  it  is  clear  that  the  only- 
significant  terms  in  the  summations  involve  J  2(1),  J  q_(1)» 
JQ(1),  J1(l)  and  J2(l).  For  n  equal  to  -2,  -1,  0,  or  2,  the 
term  is  outside  the  bandwidth  of  the  final  filter.  Then  the 
signal  output  for  SFSK  modulation  is*. 


SSFSK  =  it”  cos[Ht  +  Mt)] 


J-^l)  cos[^i)Qt  +  4p(t)] 
J^(l)  cosC^Qt  +  4cp( t ) 3 
cos[4&)0t  +  4<p(t)] 


(184) 


The  average  power  of  this  signal  is  given  by 

PQ  -  0.0033  A8  (185)- 

s 

The  output  signals  and  their  average  power  have  been 
listed  in  Table  III  for  comparison.  Note  the  significant 
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reduction  in  the  average  power  at  the  output  for  MSK  and  SFSK 
compared  to  0-QPSK.  MSK  shows  a  decrease  of  20.56  dB  from 
0-QPSK  modulation  and  SFSK  has  a  decrease  of  15.78  dB  from 
the  0-QPSK  signal.  It  can  be  argued  heuristically  that  the 
noise  power  suffers  a  comparable  loss  since  the  noise  terms 
are  subject  to  the  pulse  shaping  functions  also,  but  this 
remains  to  be  'proven. 

The  signal  to  noise  ratio  (SNR)  at  the  output  of  the 
quadrupling  circuit  will  now  be  calculated  for  the  case  of 
0-QPSK  pulse  shaping,  and  an  example  will  be  given  with  rea¬ 
sonable  numbers  assigned.  In  this  case,  the  noise  output  of 
the  quadrupler  reduces  to  the  followingi 


0-QPSK 


[-A3annc  *  A\ns  +  3A2anbnncns  -  ^Aa^n2 

+  ^Ab  n^n  -  +  ^Aa  n*}  -  kb  n^ 

2  n  c  s  4cs  2  n  c  2  ns 

*  Enc  *  !ns]  cosC1*wc't  ♦ 

+  4V0  ♦  A\ns  -  iAVcns  -  2Aanncns 


+  |A2a  b  n^  -  42a  b  n2  +  ^Ab  n^  +  ka  n^ 
2  nnc  2  nns  2nc  2ns 


+  lncns  “  |ncns  sin^ct  +  ¥*)]  (186) 


The  noise  power  is  found  by  squaring  this  and  taking  the 
expected  value  with  respect  to  the  statistics  of  nQ  and  ng. 
To  evaluate  this  expected  value,  equation  (52)  was  used  with 
W  equal  to  the  noise  bandwidth,  and  the  following  relation 
was  used  (Ref.  14i147)i 
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TABLE  III 


( 


Quadrupling  Circuit  Output 

Signal  and  Average  Power. 

Modulation 

Output 

Average 

Type 

Signal 

Power 

0-QPSK 

-A4 

2 

cos[4d)0t  +  4tp(t)] 

0.1250  A8 

MSK 

cos[4o)0t  +  4cp(t)3 

0.0011  A8 

SFSK 

w 

>  ^A4  1 

+  Ji ( 1 ) J  cos[4{i)0t  +  4<p(t)] 

0.0033  A8 

E£nm(t)} 


1-3*  •  •(m-l)(N0W)m//2 
0 


m  even 


m  odd 


(18?) 


In  addition,  it  was  assumed  that  a  ,  b  .  n  (t),  and  n  (t)  were 
mutually  independent  of  each  other  and  zero  mean.  Performing 
the  operations  gives  the  following  as  the  noise  power  output: 

256  A6(NqW)  +  720  A4(NQW)2  +  448  A2(NQW)3  +  43(NQW)4 
PN  =  I5B 

(188) 


The  signal  to  noise  power  ratio  for  the  case  of  0-QPSK  modu¬ 
lation  is  given  by 


(*) 


16  A 


8 


0-QPSK 


256  A6NqW  +  ?20  A \n0W)2+  448  A2(NQW)3+  43(NQW)4 

(189) 
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Suppose  the  intercept  receiver  is  located  10  miles 
from  the  transmitter  whose  effective  radiated  power  is  100 
watts  and  that  the  transmission  frequency  is  250  MHz.  The 
received  signal  power  is  given  by  (Ref.  13:12) 

Pr  =  PtGtGr(s£d)  (190) 

where  Pr  is  the  received  power;  P^.  is  the  transmitted  power; 
G^  and  Gr  are  the  transmitter  and  receiver  antenna  gains, 
both  assumed  equal  to  1  in  this  example;  \  is  the  wavelength, 
and  d  is  the  distance  between  transmitter  and  receiver.  Then 
the  received  signal  power  is 


Pr  =  3.5208295  x  10‘9  watts  =  -54.53  dBm  (191) 

2 

Since  the  signal  power  at  the  input  is  A  ,  then 

A  =  5.933?  x  10-5  volts  (192) 

The  noise  power  bandwidth,  W,  is  16  KHz  in  this  example  and 
(Ref.  13:13) 


(193) 


where  k  represents  Boltzmans  constant  (1.38  x  10-2^  w/Hz-°K) 
and  T°  is  the  receiver  noise  temperature  in  degrees  Kelvin 
(°K).  If  this  is  assumed  to  be  room  temperature  (290°K), 
then 

NqW  -  6.4032  x  io“17  watts  (194) 


Then  the  quadrupling  circuit  output  SNR  is,  for  this  example: 

a  3.437  x  106  =  65.36  dB  (195) 


(6), 


0-QPSK 


95 


Thus,  for  this  example,  the  SNR  is  quite  large.  If  the 
noise  power  output  for  the  case  of  MSK  and  SFSK  modulation  is 
equivalent  to  the  0-QPSK  case,  then  their  SNR's  would  be  on 
the  order  of  4 5  dB  and  50  dB  respectively,  both  of  which  are 
still  large  and  more  than  adequate  for  the  purposes  of  this 
thesis . 


Appendix  Bs  Format  A  Decoder  Calculations 
for  Two  Transmissions 

In  this  appendix,  the  outputs  of  the  Format  A  Decoder, 
shown  in  Figure  12,  will  be  calculated  for  the  case  where  two 
simultaneous  transmissions  are  received.  The  results  will  be 
evaluated  for  0-QPSK,  MSK,  and  SFSK  modulation. 

The  input  signal,  x(t),  to  the  Format  A  Decoder  is 
given  by  equation  (79).  After  mixing  with  coherent  carriers 
and  filtering,  the  signals  ares 

xj(t)  =  +  |A2an2p(t-nTc-At) 

+  §A2l>n2P(t-n'rc-Tc/2-At)  +  K(t)  (196) 

*5<t>  ■  Kv(*-”W2)  +  iA2VP(t-nTo-V2-At) 

-  |A2an2P(t-nTc-At)  -  |ns(t)  (197) 

The  signals  are  processed  thru  the  square  law  devices  to  give 

2 

ul(t)  =  ^F2(t-nTc)  +  g  A2p2(t-nTc-At) 

2 

+  |  A2p2(t-nTc-Tc/2-At) 

+  lAlA2anlan2P(t-nTc)  P(t-nTc-At) 

+  lAlA2anltn2p(t-nTc)  p( t-nTc-Tc/2-At ) 

+  ¥A2an2bn2p(t'nVAt)  P(t-nTc-Tc/2-At) 


(198) 


*  [Kanlp(t-nTc>  +  |A2an2P(t-nT0-At) 

+  §A2Vp(t-nVV2-At>]nc(t>  +  5nc(t) 

and 

2 

uQ(t)  =  jA2p2(t-nTc-Tc/2)  +  g  A|p2(t-nTc-Tc/2-At) 

2  • 

+  g  A2p2( t-nTc-At) 

+  2AlA2bnlbn2p(t-nVTc/2)  PU-nT^T^-At) 

"  2AlA2an2bnlp(t-nTc~At)  PCt-^VV2) 

'  ¥A2an2bn2p(t_nTc~At)  P(  t-nTc-Tc/2-At ) 

-  [lAlbnlp(t-nVV2)  +  iA2bn2p^-nTc-Tc/2-At) 

"  lA2an2p(t-nTc-At)]ns(t)  +  3ns(t> 

For  O-QPSK  modulation,  this  becomes 


ui(t)  =  iAi +  lAi^2  +  p2) 


+  2AlA2anlan2  +  lAlA2anlbn2 

+  ^^A2n  >) 

2  A2an2n2 


5Alanl  +  2A2an2  +  lA2an2]nc(t)  +  5nc(t) 


and 


ua(t)  =  Ja2  +  £a2(<i2  +  32) 


(199) 


(200) 
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+  2AlA2bnlbn2  '  lAlA2an2bnl 
*  ¥A2an2bn2 


► 

-  h 
2  1 

► 


b  ,  +  §A0b 


_  fi 


A0a„ 


nl  2  2  n2  "  2n2n2 
The  outputs  are  then 


:]ns(t>  *  H 


(t) 


yI(t) 


=  2AlA2anlan2  + 


[lAlA2anlbn2] 
[¥A2an2bn2]LP  +  5nc(t)] 

[(lAlanl  +  2A2an2  +  lA2an2 


)nc(t)> 


LP 


and 


yQ(t)  =  fVaW 


-  [|AiA2an2bnl]Lp 

-  *  [V^^Lr 

-  [(k 


b_,  +  ^A0b„„  -  ?A„a 


nl  22  n2  "  22n2 
For  MSK  modulation 


:)ns(t)] 


LP 


Uj(t) 


2^rrLt=At)^ 


+  2AlA2anlan2  C°S(t^)  C0S(1^^) 

+  |A1A2anibn2  C0S(^)  sin^h'^^ 

♦  ¥AkaV  -=(2^fi)  am(iiM) 


(201) 


(202) 


(203) 
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■  [k 


*nl  +  iv„2 


*  §Vn2  =in(^^)]n0( 

k  +  k(a2  +  B2)  +  k 


■("if*) 


j(t)  +  5nc(t) 

'  C  /J  ^  c 

§A|(a2  +  B2)  +  ^A2  cosfe*' 

V  r»  i 

+  sA|(a2  -  g2)  oos| 


*  5AlA2anlan2  [cos(t?  ’  tJ)  + 

+  !AlA2anlbn2  [sin(^  -  -  1 

an2b’ 


♦ft 
» 


.  J  -  -ft] 

>n2  sin(^Sl) 
fAlanl  oos(f0)+  2A2an2 


+  IpA^a, 


+  2A2bn2  sin 


C0StrJ+  tVn2  cos(f^ 
in^t-^l)]„c(t)  ♦  ln2(t) 


and 


uQ(t) 


k  sin2(fS)  *  |2A2  sin2^li^ii\ 

a  ft  /.  /  .  . 


A|  cos  | 


2/rr( 

\  'c  j 


lAlA2bnlbn2  sin 


-  -  -  ) 

2AlA2an2bnl  oos(2±E^iI)  sin(l;) 

^A2an2bn2  cos^^1)  sin^^ 


(204) 
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( 


-  [lAlbnl  sin(fJ  +  tA2b„2 

-  lA2an2  ?)]ns(t)  *  5ns(t> 

=  SAi  *  5A2(a2  +  B2)  -  hai  oos(t;) 

+  §A2(B2  -  a2)  cos  (Mj;- At  A 

+  5AlA2bnlbn2[cos^)-  cos(^  '  "T;)] 

-  IAlVn2bnl  [sin(^f;  -  X)~  Sln(^)] 

-  ^AIan2bn2  si"(2T,(^At)) 

*»©  *  lA2bn2 


2Albnl  sin' m 


-  lA2an2  oos 


(t)  +  £n2(t) 


These  results  are  then  filtered  by  the  audio  low  pass 
to  give  at  the  output  for  MSK 


yjft)  =  gA1A2anlan2  cos(%i) 

-  [lAlA2anlbn2  *ln^)]llp  + 

and 

y9(t)  =  cos (sfy 

*  [§AlA2an2bnl  ai"(^|p  *  [^"S(t)]LP 
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(205) 

filter 


(206) 


(208) 


For  SFSK  modulation,  the  analysis  is  much  more  complex. 
At  the  output  of  the  square  law  device  in  the  in-phase  chan¬ 
nel,  the  signal  is  (with  0  =  t-At)» 


Uj{t)  = 


*  2AlA2anlan2  008 


•  COS 


+  lAlA2anlbn2  008 


5*1  -  5  sin 

+  |2a|  cos2fe  -  l  sin  ^  -  ilfS] 

L  c  c  cJ 

+  J^A2  e.in2fTrt  1  4~nfl  TTAtl 

+  %  A2  sln  t"  “  sin  —  “  ~ 

*»  c  c  cJ 

+  2AlA2anlan2  co8[^  '  5  sin 

•  cosfe  -  l  sin  S|0  -  afil 

L  c  c  cJ 

+  |A1A2anl^n2  C0S[?J  "  i  sin 

1  .  4-nfl  TfAtl 

.  smlj  -  ^  sin  - 

l  c  C  CJ 

+  2§A2a  b  oos^  _  1  sin  itsQ  _  "At! 

2  A2anl Dnl  cos  T  4  sln  T  T 

L  c  c  cJ 

ai'rJfft  1  -J-  4tlO  TTAtl 

•  1  K  -  5  81n  -r;  -  —c\ 

♦  (sVnl  °°°[t0  -  5  ■*»  - 

+  IVn2  888[f  -  t  *»  ¥  -  ¥\ 

L  c  c  cJ 

♦  lA2bn2  -  5  8ln  T;  -  ^|)»o(t> 


*  2^A2anlbnl  008 


•  sin 


+  2A2bn2  8in 


+  Jn2(t) 


(208) 
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I 


Expanding  this  gives  (with  w  =  2  /T  ) 

c  c 


2  n2 

uj(t)  =  +  I  A2  +  I  A2  +  IA1  cosJwct  -  |  sin(2ooct)] 

jy  2  ^  - 

+  B  A2  cos[wc't  -  J  sin(2(1)cn)Jcos(wcAt) 

2 

+  §  sin[coct  -  |  sin(2(i3cQ)Jsin((i)cAt) 

"  i  A2  cosKt  ■  \  sin(2wcn)Jcos(wcAt) 

-  |  Ag  sin^wct  -  i  sin(2wcQ)Jsin((j)cAt) 

+  2  AlA2anlan2  cos[(wct/2)  -  £  sinUt^t)] 

•  cos[(wct/2)  -  J  sin(2(DcQ)]  cos^^^ 

2 

+  2  AlA2anlan2  cos[(wct/2)  -  £  sin(2wct)] 

.  sin[(dJct/2)  -  ^  sin(2wc0)]  sin^r^ 

+  §AiA2anlbn2  cos[^ct/2)  -  ^  sin(2d)ct)] 

•  sin[((oct/2)  -  ^  sin(2(0cn)]cos^^^ 

"  fAlA2anlbn2  cos[^ct/2)  -  \  sinUc^tj] 

•  COS  [(oct/2)  -  J  sin(2o)cQ)]  sin^y^ 

+  lFA2anlbn2  j^cos(^ct/2)  -  \  sin^O)]  cos(t^) 
+  sin[((»)ct/2)  -  J  sin(2ti)cfl)]  sin^^^j 


2  A2anlDn2 


+  sin 


« 


I 


-  cos 


j^sin[((oct/2)  -  J  sin(2wcQ)j  cos^^ 

[(wct/2)  -  J  sin(2(,)cn)]  sin^^J 

lAlanl  cos[^ct/2)  -  J  sin(2(oct)] 

+  2A2an2  cos[(wct/2)  -  £  8in(2u)cfl)]  cos^^ 

+  2A2an2  sin[(wcV2)  -  J  sin(2wcn)] 

+  fA2bn2  sin[(wct/2)  -  £  sin(2cocn)j  cos(^j 


"  fA2bn2  cos[(wct/2)  -  J  sin(2wcQ)]  sin^ 


nc(t) 


+  K(t) 


(209) 

The  sinusoids  are  now  expressed  as  Bessel  function  expan¬ 


sions. 


Uj(t) 


2  2 

1A1  +  I  A2  +  8  A2  +  IA1  f  J i ( i/2 )  cos(l+2i)(l)ct 

2  2 

+  8  A2  cos(^ca^)  ?  J ( 1/2 )  cos(o)ct  +  2i(i)cQ) 

2 

+  i  A2  sin(^cAt)  ?  J ( 1/2 )  sin(d)ct  +  2i(ocfl) 

"  f  A2  C0S(f,JcA't)  £  ( 1/2 )  cos((oct  +  2i(i)cQ) 

q2  £ 

“  8  A2  sin(wcAt)  I  J± ( 1/2 )  sin(wct  +  2i(ocn) 

+  003 (t^)  |  JiUA) 

•  2  Jj(l/4)  cos[(wct/2)  +  2jwcn] 
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T 


+  2AlA2anlan2  Sin^T^)  i  GOa(^jrWGty 

•  E  J.(l/4)  sin[(wct/2)  +  2jc*)c0j 

+  2AlA2anlbn2  C0S^^)i  cos^Tici)ct) 

•  2  J^d/4)  sin[((i)ct/2)  +  2jwc0l 
«j 

-  2AlA2anlbn2  Sin(~T^)  i  cos(^T^^c^ 

•  E  Jj(l/4)  cos  [(wct/2)  +  2jwcQ] 

J 

+  ¥A2anl  bn2  C°s2(D^)  l  Ji(l/4)  cos[(a>ct/2)  +  2iwcn] 

•  E  Jj(l/4)  sin[((i)ct/2)  +  2jwcn] 

‘  fA2anlbn2  cos(^)  sin(^) 

•  E  J.^l/4)  cos[(wct/2)  +  2i(Dcn] 

•  E  J. (1/4)  cos  [(d)ct/2)  +  2jwca] 

J 

•  ¥A2anlbn2  cos(t;)  sin(^) 

•  E  J^l/4)  sin[((i)ct/2)  +  2i(i)co] 

•  E  J^d/4)  sin[(03ct/2)  +  2jwcn] 

J 

-  ^A2anibn2  sin2^-p^  ?  J^l/4)  sin  [(coct/2)  +  2icoc0] 

•  E  Jj(l/4)  cos[(d)ct/2)  +  2jd>c0] 

4  Kara  I  Ji<1/'*>  00s(iri<"ct) 

+  2A2an2  C°S(^)  f  Ji^A)  cos[(d)ct/2)  +  2i(»>cn] 
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+  lA2an2  sin(^jl  Ji^A)  sin[(uct/2)  +  2iwcn] 

+  lA2bn2  cos^T^)i  J^dA)  sin[(d)ct/2)  +  2iwcn] 

-  |A2bn2  Sln(^)  i  Ji(1A)  cos[(wct/2)  +  2iwca]jnc( t) 

+  ^(t)  (210) 

A  careful  examination  of  equation  (210)  shows  that  very  few  of 

these  terms  pass  the  output  filter.  The  filtered  output  is: 

yI(t)  =  5AlA2anlan2  C°S(~T^)  Jq^A) 

"  [lAlfl2an2bnl  sln(“T~)  Jo(  1,/4']Lp 

+  [K^LP  l211> 

Through  an  identical  analysis,  the  quadrature  channel  output 
is  found  to  be: 

ya(t)  =  gA1A2bnlbn2  cos^)  J2(l/4) 

+  [kV^nl  J0(lA)]Lp 

♦  [K(t,]lP  <212> 
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